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We investigate the influence of the thermal properties of the boundaries in turbulent 
Rayleigh-Benard convection on analytical bounds on convective heat transport. Using 
the Doering-Constantin background flow method, we systematically formulate a bound- 
ing principle on the Nusselt-Rayleigh number relationship for general mixed thermal 
boundary conditions of constant Biot number rj which continuously interpolates between 
the previously studied fixed temperature (77 = 0) and fixed flux (rj = 00) cases, and 
derive explicit asymptotic and rigorous bounds. Introducing a control parameter R as 
a measure of the driving which is in general different from the usual Rayleigh number 
Ra, we find that for each 77 > 0, as i? increases the bound on the Nusselt number Nu 
approaches that for the fixed flux problem. Specifically, for < < 00 and for suf- 
ficiently large R {R > Rs — 0{rj^^) for small 77) the Nusselt number is bounded as 
Nu < c{ri)R^/^ < CRa^'^, where C is an 77-independent constant. In the R ^ 00 limit, 
the usual fixed temperature assumption is thus a singular limit of this general bounding 
problem. 



1. Introduction 

Rayleigh-Benard convection, in which a fluid layer sandwiched between two horizontal 
plates is heated from below, has long attracted considerable experimental and theoretical 
attention; this is due not only to its importance as a model with numerous applications 
in engineering, geophysics, astrophysics and elsewhere, but also because it has proved 
such fertile ground for explorations concerning stability an d dynamics, pattern forma- 
tion and - Tinder sufficient heating — convective turbulence ([Cross fc Hohenberd (Il993l) : 



[KadanoflE (2001)). Particular attention has been given to the Nusselt number Nu, repre- 
senting the convective enhancement of vertical heat transport, and its dependence on the 
Rayleigh number Ra, a measure of the driving via the temperature difference across the 
fluid layer. This dependence appears to take a scaling form Nu ~ f{Pr, T)Ra^ (with pos- 
sible logarithmic corrections), where F represents geometric effects and Pr is the Prandtl 
number, and a major goal of theory and experiment is to find p. 

Observations that the heat transfer is essentially confined to thermal boundary lay- 
ers near the plat es, separated b y an isothermal core, suggest that p = 1/3, as pro- 
posed already by 'Malku d (Il954 ); th is prediction appears to be consistent with large- 



Ra experiments (Nicmel a fc Sreenivasa n (20066)) and numerical simulations of turbu 



lent Boussin esq convection ( Amati et al\ ([20051)). In some experiments ( Chavanne et all 



( 19971 [200l[ )) an increase in the scaling exponent was observed at the highest accessible 



Ra values, suggesting a possible transition to a p — 1/2 asymptotic regime predicted 
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bv IKraichnan ( 1962), but other experi n ients at cq n ipara bl e Ra have failed to o bserve 
such a transition ( Glazier et al. ( 19991 ): Sommerial (jl999l ): Niemela et al. ( 2000^ )). and 
possible Prandtl number variabilities or non-Boussin esq effects may have played a role 
( Glazier et "aLl (ll999l) :l iNiemela fc SreenivasanI ([2006313)); strong evidence of this so-called 
"ultimate" regime appears to be absent. 

Analytical bounds on heat transport: 

While good agreement with experimentally observed scaling for a range of Rayleigh and 
Prand tl numbers has been attained by a phcnomenological theory (jGrossmann fc Lohse 
(|200lh 1. we focus on results derived systematically from the underlying governing differ- 
ential equations. Uppe r bounding pri nciple s derived und er some statistical assumptions, 
dating to the work of iHowardI (1 19631) andlBussd (Il969l), have yie lded Kraichan's expo- 
nent p = 1/2. More recentlv. iDoering fc ConstantinI ( 1992 . 19961 ) realized that an idea 
of decom posing flow variables into "background" and "fluctuating" components, intro- 
duced bv lHopj ( 194lh to prove energy boundedness, could be extended to obtain rigorous 
analytical bounds (without additional assumptions) on bulk transport quantities. The 
Doering-Constantin "background flow" method has since proved remarkably fruitful in 
obtaining bounds in good quantitative agreement with experiment or direct numerical 
simulation in a wide range of flows. 

For Rayleigh-Bcnard convection with fixed temperatures at the lower and upper bound- 

1 /'? 

aries of the fluid, the backgrou nd flow method yields a rigor ous bound Nu < Cq Ra ' 
uniform in Prandtl number Pr ('Doering & Constantin and while extensive sub- 

sequent investigations (,Ker swcll (1997, 2001); Plasting & KcrswcU (2003)) have improved 
and optimized the constant Cq in the bound, for general Pr it has to date only proved pos- 
sible to lower the exponent p fror n the Kraichnan value p = 1/2 under add i tional length 
scale or regularity assumptions (jConstantin fc Doeringl (jl99 B); iKerswelll (|200lh . The 
assumption of infinite Prandtl number, however, by imposing an additional constraint 
on the velocity field, permits a lower i ng of t he scaling exponent to jo = 1 / 3 (with pos- 
sible logarithmic corrections) ( Chan ( 197lf ): Constantin fc Doering (| 19991 ): Her lev et al. 



2006h'): the best curr ent rigorous bound in this case has the form Nu < C^lnRaW^Ra}^ 



Doering et al. ( 20061 )). and related results have recently been obtained bv IWang ( 20081 ) 



for sufficiently large finite Prandtl number. 

Influence of thermal properties of the plates: 

The above analyses were performed under the usual assumption that the lower and 
upper boundaries of the fluid in Rayleigh-Benard convection are held at known uniform 
temperature, or equivalently, that the bounding plates are perfect conductors. In prac- 
tise, though, the boundaries are imperfectly conducting; and the thermal properties of 
the boundaries have long been understood to affect the initial instabilit y to convection 
and the weakly n onlinear behavio ur beyond transition ( see for in stance Sparr ow et al. 
( 1964) : iHurle et al. (1967): C hap man fc Proctol (|l980l ): iBusse fc R iahi (1980) and the 



review m 



Cross fc HohenbCTi (il993h ). Even when the bounding plates have much higher 



conductivity than the fluid, as the Rayleigh number (and hence the Nusselt number) 
increases, the effective conductivity of the fluid, depending on Nu, eventually becomes 
comparable to and then exceeds that of the plates; in fact, in the Ra oo limit one 
might expect the fluid effectively to "short circuit" the system, with the bounding plates 
acting essentially as perfect insulators by comparison. 

The effect of the finite thermal conductiv ity of the bound i ng pla t es on convective he at 
transport has stimulated recent m odelling IChauniat et a/.l ( 2002 ): [Chilla et al. ( 2004f )). 
experimental ( Brown et al. ( 20051 )) and numerical (jVerzicco ( 2004 )) studies with the aim 
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of reconciling var ious expe rimental results with each o ther and with fixed -temperature 
theoretical ( Grossmann fc Lohsc (j200i)) and numerical ( Amati et al. ( 20051 )) predictions. 
Recent numerical simulations comparing fixed fiux and fixed tem perature boundary con- 
ditions have reached differing conclusions: The computations of IVerzicco fc Sreenivasan 
([2OO8,) in cylindrical geometry found Nu ~ Ra}^^ scaling, but that for a given large 
enough Ra, the Nusselt number is reduced upon replacement of fi xed temperature con- 
ditions at th e lower boundary of the fluid by fixed flux conditions. iJohnston fc Doering 
(|2007l . 120081 ). on the other hand, found that in numerical integration of two-dimensional, 
horizontally periodic convection, the heat transport for large Ra was the same, namely 
Nu ~ Ra'^^'' , for fixed temperature and fixed flux conditions at the upper and lower 
boundaries of the fluid. 



Pre dating most of the above recent investigations and with similar motivations. I Otero et 
initiated the analytical study of the effects of thermal boundary conditions on 
Rayleigh-Benard convection using the background flow method, obtaining a bound on 
the heat transport with fixed flux boundary conditions at the fluid boundaries, which 
again took the for m Nu < CooRa^^^ (this work was recently extended to porous medium 
convection by Weil (2007)). Howev er, the mathematical structure of the flxed flux bound- 



ing problem of lOtero et 



( 2002h . and various intermed iate scaling results, tur n ed ou t 
to be qu i te diff erent from the fixed temperature case ( Doering fc Constantin ( 19961 ): 
Kerswell (|200ll) ). It is thus natural to wonder how these two extreme cases, correspond- 
ing respectively to the idealizations of perfectly conducting and insulating plates, are 
related vis-d-vis their bounding problems, and which is more relevant to real, finitely 
conducting boundaries. 

Outline of this paper: 

In the present work we reconsider the effect of general thermal boundary conditions on 
systematica lly derived analytic al bounds on thermal convection, continuing the program 
initiated bv lOtero et al. (2002); for simplicity we consider only identical thermal proper- 
ties at the top and bottom fluid boundaries in the mathematically idealized horizontally 
periodic case. We consider a common model for poorly conducting plates, namely mixed 
(Robin) thermal boundary conditions of "Newton's Law of Heating" type, with a fixed 
Biot number 77, so that ry = gives the fixed temperature and rj = 00 the fixed fiux case; 
to our knowledge th e only prior bounding study with general Biot number is the work of 
Siggers et al\ (|2004l ) on horizontal convection, in which mixed thermal conditions were 
imposed at the lower boundary of the fiuid. 

In Section [2] of this paper, we carefully develop a general formulation and bounding 
principle using the Doering-Constantin background flow method, and in Section [3] spe- 
cialize to mixed thermal conditions in a manner that interpolates smoothly between the 
fixed temperature and fixed flux cases. The use of a piecewise linear background tem- 
perature profile and exphcit estimates derived in Section [4] enables us, in Section [U to 
derive analytical bounds on the Nu-Ra relationship asymptotically valid for Ra —f 00. 
For completeness, we also prove the boundedness of temperature and velocity fields in 
Appendix |^ and prove rigorous, though less sharp, bounds on the Nu-Ra relationship 
in Appendix [Cl 

Summarizing our results: Since in general the boundary temperatures are unknown 
a priori, it is necessary to introduce a control parameter R, which equal s the standard 
Rayleigh number Ra only in the fixed temperature case r] = O. IOtero et al ] 12002) showed 
in the fixed flux case = 00 that while the bound was Nu < CooRa^' , it was obtained 
through the estimates Nu < ciR^/^, Ra > C2R^^^, quite unlike the fixed temperature 
case Nu < Coi^^/^ Ra = R. 
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In the present work, for general sufficiently small Biot number 77 we show that for 
small R we have Nu < 0{R^/^), Ra > 0{R) as in the fixed temperature case, but for R 
(and hence Ra) beyond some critical parameter Rg — 0{r]^^), we find Nu < ci{jf)R^/^, 
Ra > C2{ji)R^/^, implying Nu < C^iRa^^"^ with intermediate scaling as in the fixed flux 
case. Interestingly, for 77 > we find = Coo- at least at the level of our estimates, 
the asymptotic scaling in each case is as for fixed flux boundary conditions (providing 
rigorous support for the intuition that for sufficiently high iZa, the plates essentially act 
as insulators), while the fixed temperature problem is a singular limit of the general 
asymptotic bounding problem. More details of the scaling of the bounds in different 
regimes, together with numerically o btained conservative bo unds for piecewise linear 
backgrounds, will be given elsewhere ([Wittenberg fc "Gaol EOOH )). 

The use of mixed "Newton's Law of Cooling" boundary conditions with fixed Biot num- 
ber r] to model imperfectly conducting boundaries is a simplification, how ever, since in 

genera l the Biot number depends on horizont a l wave number (see for instance lCross fc Hohenbere 

] n a subsequent paper ( IWittenberd (|2008f) . l^art 2 of the present work), we im- 
prove upon our model by formulating and obtaining bounds for the more realistic problem 
of a fiuid bounded by plates of finite thickness and conductivity, establishing a systematic 
correspondence between that situation and the present fixed Biot number case. 



2. Governing equations and bounding principle with general thermal 
conditions at fluid boundaries 

We begin by formulating the standard Rayleigh-Benard convection problem in the fluid 
and developing a bounding principle in the usual way, but without fixing the thermal 
conditions at the fiuid boundaries. For reference and clarity, though, we occasionally 
point out the forms of our results in the fixed temperature and fixed fiux special cases 
previously treated in the literature, as our development is designed to interpolate between 
these extremes. 

2.1. Governing differential equations and nondimensionalization 
In the Boussinesq approximation, the equations of motion in the fluid are 

^ + u* • V*u* + —V*P* = UfW^n* + agiT* - To)e^ , (2.1) 
at* Pf 

V* • u* = , (2.2) 
u* • V*T* = KfV*^T* , (2.3) 
u*U,^o,h = 0, (2.4) 



dT* 



where Vf and Kf are the momentum and thermal diffusivities of the fiuid, respectively, g is 
the acceleration due to gravity, a is the thermal expansion coefficient, pf the fluid density 
at some reference temperature Tq, and h is the height of the fluid layer, as in Figure [TJ 
we also let Cpj be the specific heat and — Pf Cpj Kf he the thermal conductivity of 
the fluid. In this formulation, the compressibility of the fluid is neglected everywhere 
except in the buoyancy force term, and the pressure is determined via the divergence- 
free condition on u* . Variables with an asterisk are dimensional, we have no-slip velocity 
boundary conditions in the vertical direction, and we take periodic boundary conditions 
in the horizontal directions, with periods L* and L*, respectively. 

The present formulation is designed to be flexible with respect to the choice of thermal 
boundary conditions (BCs) at the plate-fluid interfaces z* = and z* = h; so in general 
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z' = h 



z- = 



Upper thermal BCs 




Fluid 



Lower thermal BCs 



Figure 1 . Geometry of Rayleigh-Benard convection system with thermal boundary conditions 
imposed at upper and lower limits of fluid layer. 



we let be a given temperature scale, determined according to the thermal BCs, and 
introduce a reference ( "zero" ) temperature Trof • 

The usual and most-studied assumption regarding thermal boundary conditions at the 
interfaces is that the temperature is fixed at the upper and lower boundaries: 

T%.^^ = T*{x\y\QX)^T;, r*|,.^,, = t;. (2.5) 

These Dirichlet BCs imply a natural choice of reference temperature Trof = Jt*, while the 
imposed temperature drop AT* = — T*\^,^q = T^—T^ introduces a natural temperature 
scale Q = AT* . At the opposite extreme is the fixed flux assumption that the thermal 
heat flux —XfT*, = —A/ dT* /dz* through the fluid boundaries is a constant, which we 
call This corresponds to the Neumann BCs of fixed normal temperature gradient —(3* 
at the interfaces: 

t;. |,-=o = t;. = -/?* = (2.6) 

the corresponding temperature scale is = h(3* = h^/Xf, while in this case Trcf is 
arbitrary. More general mixed (Robin) thermal BCs are discussed later in Section [3l 

We nondimensionalize using Trcf and the temperature scale Q, and with respect to the 
fluid layer thickness h and thermal diffusivity time j kj; that is, we take h, h'^/Kf, U = 
Kf/h and PfU^ as our appropriate length, time, velocity and pressure scales respectively. 
For Trcf 7^ Tq, the nondimensional fluid momentum equation will contain a constant 
term proportional to Trof — Tq hi the direction, which we can take care of by absorbing 
it into the rescaled pressure; this effect of the buoyancy force corresponds to a linear 
vertical pressure gradient. In summary, the nondimensional variables (without asterisks) 
are deflned by: 

X* ^ t* u* ^ T*-Tref 1 P „Trcf-Tn 

h tsced U e Pr pfU^ 6 

where tsced — h'^/i^f, U — Kf/h, and Pr and R are defined below. The dimcnsionless 
periodicity lengths in the transverse directions are — L%/h and Ly — L*/h. 

The equations for the nondimensional fluid velocity u = (u, v, w) and fluid temperature 
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T are thus 



Pr-^ + u-Vu] +Vp = V^u + RTe^, (2.7) 



V-u = 0, (2.8) 
BT 

u- Vr = V^T, (2.9) 



dt 

with no-shp BCs uj^^g,! = Oj L^, Ly-periodic BCs in the horizontal x and y direc- 
tions in all variables. Here the usual Prandtl number Pr and the control parameter R 
are defined by 

Pr = ^, R=^e. (2.10) 

Kf VfKf 

In the fixed temperature (Dirichlet) and fixed flux (Neumann) limits, the nondimensional 
thermal BCs are, respectively, 

r|.=o = l, 7^1.=! =0, (2.11) 

Some additional notation: 

Following Otero et all ( 2002f ). for functions h{x, y, z) and g{t) we define the horizontal 
and time averages, h{z) and {g) respectively, by 



and 



h{z) = ^ j j h{x, y, z) dxdy = ^ ^ j K^, V, z) dxdy 
{g) = limsup- / g{t)dt, 



where A — L^Ly is the nondimensional area of the plates. Also, J h ~ A h,{z) dz 
denotes a volume integral over the entire fluid layer, and the norms arc defined over 
the fluid by 

= J ^ J jj {x , y , z) dxdy dz . 

Finally, over the domain we consider with horizontally periodic BCs, surface integrals 
of vector fields H over the fluid boundary (the interfaces T between the plates and the 
fluid) are 

1 1 



nUdS ^ AUe^ 



2.2. Rayleigh and Nusselt numbers 

We define the nondimensional horizontally- and time-averaged temperature drop across 
the fluid as 

= -{T\\j = {TU=o TU=i) = (2.13) 

where AT* — {T*\z*=o — T*\z'=h)] we observe that this is known a priori only for fixed 
temperature BCs, in which case AT* = T^ — T^ — 6, AT — 1. The conventional Rayleigh 
number Ra is defined in terms of this averaged temperature difference as 

Ra=^AT*^^^AT, (2.14) 
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showing that Ra is related to the control parameter R by 

Ra^RAT. (2.15) 

The Nusselt number Nu is a nondimensional measure of the total heat transport 
through the fluid layer, which may be derived by first rewriting the thermal advec- 
tion equation in the fluid (|2.9p as a conservation law, + V • J ~ (using incom- 
pressibility). Here the dimensionless heat current J = uT + J^, is composed of the 
conductive heat current — — VT and the convective heat current J„ = uT, with 
corresponding overall instantaneous conductive and convective vertical heat transport 
jc{t) = /ez • Jc = /X4 Jq ~Tz dzdxdy = -A T\_^^^ and jy{t) = /e^ • J^, = / wT, respec- 
tively. The Nusselt number Nu is now defined as the ratio of the total (averaged) vertical 
heat transport, j{t) = J ■ J = jc{t) + jv{t), to the purely conductive transport: 

A more useful expression, which allows us to estimate Nu from the equations of motion, 
is found by relating (J wT) to the time-averaged temperature drop and boundary flux. 
To do so, we begin by taking the horizontal average of the temperature equation (|2.9p . 
using the horizontally periodic BCs, to get 

rir Fir f) 

Integrating over z and using the vertical no-slip boundary conditions on w, 

ijT^Aj\^- T.)^ dz^ljT + A (-n)l t„ = 0. (2.18) 

As we show in Appendix \X[ the thermal energy ||r|| = JT^ is uniformly bounded in 
time for the thermal BCs we consider, so that J T is also uniformly bounded by the 
Cauchy-Schwarz lemma via J T < A^l"^ ||T||. Hence on taking a time average of (|2.18p . 
the time derivative term vanishes, and we find that 

(-7^^)12=0= (-^^)L=i^/3, (2.19) 

where the above expression defines /3, the horizontally- and time-averaged vertical tem- 
perature gradient, or equivalently, the nondimensional heat flux at the fluid boundaries; 
note that (only) in the fixed flux case, this quantity is known, f3 = 1. As expressed by 
(|2.19p . obviously on average, there is a balance between the heat fluxes entering the fluid 
layer at the bottom and leaving at the top. 

Taking a time average of (j2.17p . via a maximum principle on T the time derivative 
term (dT/dt) would vanish; ho wever, for fix e d flux BCs we do not have such an a priori 
maximum principle. Following Otero et al. ( 2002h . uniformly in thermal BCs we can 



instead multiply (|2.17p by z and integrate to obtain 



d_ 

di 



1 _ /•! _ 

zTdz+ z{wT -Tz)^ dz = ; (2.20) 



/o Jo 

integrating the second term by parts and using the no-slip boundary conditions, we have 

J z {wT - Tz) dz = z {wT - Tz)\l- J {wT - T^) dz = -Tz\z=i--^ J wT+T\1^^ . 

(2.21) 
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Now as before, via 

and the uniform boundedness of ||r||, the time average of the first term in (I2.20p vanishes. 
Thus taking time averages of (|2.20|) . using (|2.2ip and the definitions (|2.13p and p.l9p . 
we obtain 

^(^Jwt'^^/3~AT. (2.22) 

The expression (|2.22l) may now immediately be substituted into (|2.16p to obtain the 
fundamental identity for the Nusselt number for general thermal BCs, 

^« = (2.23) 
and consequently, using (|2.15p . we have Nu Ra ~ Rj3. 



In t he special case of fixed temperature BCs, the identity (|2.23p is well-known (jPoering k, Constantin 



( 1996( )): with AT = 1, the Nusselt number equals the time- and horizontally- averaged 



flux at the boundary, Nu = (3 = {—T z\z=o.i) i thus an upper bound on Nu is found 
by bounding [3 from above. For the fixed flux case, /3 = 1 is k nown, the Nu s selt n umber 
is the inverse of the averaged temperature drop, Nu — AT~^ ( Otero et al. ( 2002[ )). and 



one seeks a lower bound on AT. In general, though, neither AT nor (3 is known a priori, 
though they are related through the thermal BCs, as discussed in Section [3l 

2.3. Global energy balance 

We next obtain the basic "energy" identities from the governing Boussinesq equations, 
which allow us to relate the Nusselt number to the momentum and heat dissipation: 
First, taking the inner product of the momentum equation (|2.7I) with u, integrating over 
the fluid domain, integrating by parts and using incompressibility and no-slip BCs, we 
find 

^ "uf = - II Vu||' + R fwT. (2.24) 



2Pr dt 

The norm of the fluid velocity is a priori bounded, as shown in Appendix [Xl hence, 
taking time averages, we derive the result (using (|2.22p ) 



^{\\Vuf) = QwT'j=AiP-AT). 



(2.25) 



Observe that (I2.25P implies that /3 > AT, so that by p.23p we have Nu > 1, as expected. 

Similarly, we can take multiply the thermal advection-diffusion equation (|2.9p by T, 
integrate over the fluid and integrate by parts. Neither the advection term nor the hori- 
zontal boundary terms contribute, so we flnd 



J n • (TVT) dS = - \\VT\f + A TT,\l^^ . (2.26) 
Taking time averages and using boundedness of ||T||, we find 

(IIVTII^) =a(t7;|;^o). (2.27) 
Again, for reference we quote these results in the known limits : for fixed temperature 



BCs (|2.1ip the global energy identities (|2.25p . (|2.27p become (jPoering fc Constantin 



Bounds for convection with fixed Biot number boundaries 

(1991)) 



-l^{\\Vnf) = {fi-l) = Nu-l, (2.28) 
l(|lVr||2) = (-7r(0))=/3 = iVw; (2.29) 



while in the fixed flux case (|2.12|) . they are ( Otero et al. ( 20021 ) ) 



l-{\\Vu\\^)^{\~m = l-Nu-\ (2.30) 

VT||^) ^ {-T\\j - AT = Nu-\ (2.31) 

2.4. Background fields 

In the spirit of the "background" method of Doering and Constantin, we now introduce 
a decomposition of the flow into a background, which carries th e boundary co nditions 
of the flow, and a spac e- and time-dependent fluctuating fleld ( Doering fc Con stantinI 



(1991); iKersweUl (|2n0lh ). One chooses flelds U and r which satisfy the same velocity 



and temperature BCs as u and T, and appropriate evolution equations — such (U, t) are 
referred to as "background" flow and temperature fields — and lets v(x, t) and 0(x, t) be 
arbitrary space- and time-dependent perturbations satisfying homogeneous BCs, so that 
the velocity and temperature fields are decomposed into a background plus a fluctuation, 
according to u(x, t) + v(x, i), T(x, t) — t + 9{x, t). 

For simplicity of presentation we restrict our attention to the case of zero background 
velocity field and a z-dependent temperature background, (U, r) = (0,r(z)), as it ap- 
pears th at more general backgrounds are unlikely to improve the overall scaling of the 
bounds (iKerswelll (I2OOII )). Furthermore, when the upper and lower boundaries of the 



fluid have identical thermal properties, it is sufficient to consider only background flelds 
satisfying r'(0) = t'{1) (compare (|2.19p l. and we deflnc 

Ar = r(0) -t(1), 7 = -r'(0) = -t'(1). (2.32) 

We thus define v and 6 via the decomposition 

u(x,t) = v(x,i), T(x,t) = r(z) + 0(x,i), (2.33) 

(note that we prefer to preserve the (notational) distinction between the overall velocity 
field u and the fluctuating fleld v, though for convenience we ignore this distinction in 
writing the components (u, w, w) of the velocity field, for instance w = u • e^. = v • e^) and 
immediately obtain identities between the norms of the gradients of the full solutions 
and their fluctuations: 

l||Vuf = l|lVvf , (2.34) 



II VT|r = II Veir + 2j e,y + j t'\ (2.35) 

Inserting the decomposition (j2.33p into the Boussinesq equations yields the evolution 
equations for the perturbations, 

Pr"^ ^l^+v- Vv^ +Vp = V2v + i?6'e^, (2.36) 

V-v = 0, (2.37) 

-X- 1 

dt 



f)0 

— +v V6l = V26l + T"-wr', (2.38) 
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where v and satisfy appropriate homogeneous boundary conditions; that is, all fields are 
horizontally periodic, v satisfies the no-slip BCs v|z=o,i — 0, and 6 satisfies homogeneous 
thermal BCs consistent with those for T. In particular, for fixed temperature BCs, we 
have t(0) — 1, r(l) = 0, and 9 satisfies homogeneous Dirichlet BCs, O\z=o.i = 0; while 
in the fixed fiux case, we require r'(0) — t'(1) = —7 = —1, so that 9 satisfies the 
homogeneous Neumann BCs 6*^12=0,1 = 0- 

In an analogous way to the calculations of Sect ion [2. 3 1 we may find the energy identities 
for the fiuctuations v and 9: Taking the inner product of (|2.36p with v, integrating over 
the fluid and using incompressibility, the evolution equation for the norm of the 
fluctuating velocity field is 



1 



2PrR dt 



1 

R 



w9, 



(2.39) 



as in (|2.24p . We find the corresponding evolution for the perturbed temperature 9 by 
multiplying (|2.38p by 9 and integrating: 



v4 



lz=0 



^j9t'. 



(2.40) 



For future reference, using (|2.32p and the decomposition (|2.33p . the second boimdary 
term in (|2.40p has time average 



= -1{ 



= 7A6I = 7 (AT - At) 



(2.41) 



2.5. Governing equations for general bounding principle 

In order to formulate upper bounding principles for the Nusselt number, we take appro- 
priate linear combinations of the identities (|2.39p - (|2.40p for the fiuctuating quantities, 
and the identities (|2.34p - (|2.35p for the decomposition of the gradient into background 
and fiuctuating parts. In general, such linear combinations may contain three free param- 
eters, over which one might optimize to obtain the best possible bound available within 
such a formalism ( Kerswell ( 1997 . 200 1[) ). Ho wever, we shall consider only the restricted 
special case of a single "balance parameter" b (jNicodemus et al. in this simplified 

formulation, in fact we do not require the energy dissipation equation (j2.39p at all. 

We can eliminate the / 9zt' term in the thermal energy equations by taking 2 • (I2.40p -I- 
((05)) to give 



ld_ 

2di 



{2\\9f) 



IVTII 



2 / 9wt' + 2A\ 



2 A 9t'\ 



(2.42) 



Now taking a weighted average b ■ ([2^ + {1 - b) ■ (|2:34| = b ■ [2 (|2^ -t- ([235| ] + (1 
b) ■ [0 (^:^1) + (11311)], we have 



2 dt 



(2b\\9f^ +b\\VT\\ 



Vl|Vu|r^5 



2b A 9t'\ 



z=0 



-2bA 99,^ 



1 - b 



R 



Vvp -2bT'w9 -b\V9\^ 



(2.43) 



We now take time averages and note that the time derivative term vanishes due to the 
boundedness of He'll, as shown in Appendix El Using ([2?25| . ([2?27l) and ^(2A\\ . we find 



&QrA[v,0], 

(2.44) 



bA{TT,\^^^) + b)A{p - AT) = 6 / t'^ + 2bA'yA9 + 2bA{99,\^^^ 



Bounds for convection with fixed Biot number boundaries 
where we define the quadratic form 



2-iiJv,e] = (^j 



1, 



hR 



IVvl- 



2t'w0 



1 



Here we have defined an "effective control parameter" via 



Rr — 



1 
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(2.45) 
(2.46) 



having observed that the quadratic form Qt.Ro depends on R and the balance parameter 
b only through the combination bR/{b— 1). We desire a positive balance parameter b so 
that a lower bound on Qr.Ra (and hence on Q'^ ^ in (|2.48p - (|2.49p below) should imply 
an upper bound on /3 and/or a lower bound on AT; since a necessary condition for QT,Rg 
to be a positive definite quadratic form is that i?o > or (6 — l)/b > 0, we thus require 
6> 1. 

Continuing with the formulation of the governing equations, using (|2.33p . (|2.19p and 
p.32p . we decompose the first term in (j2.44|) via 



(ttJ 



z=0' 



z=0' 



r > 

\z=0' 



r ) 

\z=0' 



13 At + 7A6' - 



IL»> 



(2.47) 



Substituting (|2.47p into (|2.44p . writing A6 = AT — Ar and rearranging terms, we obtain 

b 



(1 - b){(3 -AT) + b {(3At - -fAT) = b (^J^ t'^ dz - 7A1 
where we have now defined the quadratic form 



A 



2t'w6 + \V6\ 



A ee. 



(2.49) 



Equation (|2.48p . which is still fully independent of thermal BCs (subject to the symmetry 
condition t'(0) = ''"'(1) — —7), is the governing identity underlying our upper bounding 
principle. 

We comment that the prime in the notation Q'^ refers to the addition of the bound- 
ary terms to Qt,r^ (no implied differentiation), and note that for both fixed temperature 

and fixed flux BCs, the boundary term (^^z|^_o) vanishes, so — Qt,r^- In these 

cases, (|2.48p thus reduces to 

»i 



Nu-l=P-l = b 



(2.50) 



for Dirichlet thermal BCs (for which AT — At = 1), and to 

1- Nu^^ = 1- AT = b(^J^ T'^dz-2AT+1^ - ^QrMA^,0]- (2.51) 

for Neumann thermal BCs, in which case /3 — j — 1. For general thermal BCs, f3 and 
AT are both a priori unknown, but they are related via the boundary conditions (see 
for instance (|3.5p below), so that (|2.48p may be written in terms only of either (3 or AT. 

2.6. Allowed fields, admissible backgrounds and the spectral constraint 

As formulated thus far, the general governing equation p.48|) is an identity. If, for given 
thermal BCs, one had access to v(x, i) and 9{x,t) satisfying (|2.36P " (|2.38p . or sufficient 
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information about them to compute the time-averag ed quadratic form Q^.^Jv, 9] ^M^, 
then f3 and AT, and hence Nu, could in principle be computed. However, for turbulent 
convection such analytical information is well beyond the limits of what is (currently) 
accessible. The fundamental insight underlying upper bounding methods for convection is 
firstly, that if for some r and b, Q'^ ^ can be shown to be bounded below, then this yields, 
ultimately, an upper bound on the Nusselt number Nu (for a given R) ; and secondly, that 
such a lower bound on ^ may indeed often be demonstrated provided one is prepared 
to widen the class of fields v, over which the minimization takes place, as long as this 
class contains all solutions of (|2.36p - (|2.38p . The cost of weakening the constraints on v 
and is that this may reduce the lower bound on Q'^ ^ and thereby weaken the upper 
bound estimate for Nu. 



Allowed fields v, 6: 

In considering the class of allowable flows v = u and (fluctuation) temperature fields 
9 = T — T over which to minimize Q'^ , we observe that if the dynamical constraints 
(|2.36p - (|2.38p on the fields [v, 0] are removed without being replaced by assumptions 
on the temporal structure or correlations of the fields — detailed knowledge of which is 
unavailable — it becomes sufficient to minimize the quadratic form Q'^ ^ over stationary 
fields v(x) and 0(x). The conditions we can assume these fields v(x) and 6'(x) to satisfy 
are: appropriate homogeneous boundary conditions, the incompressibility constraint on 
the velocity fluctuations, and boundedness of ||v|| and \\9\\ (see Appendix More 
precisely, we denote the "allowed" fields v and 9 over which we minimize Q'^ ^ to be 
defined on [0, L^] x [0, Ly] x [0, 1], periodic in the horizontal directions, so that V • v = 0, 
V = on z = 0, 1, and 9 satisfies homogeneous thermal boundary conditions consistent 
with those of T . For fixed temperature convection, plausible but non-rigorous regularity 
assumptions restrictin g the class of allowed fie l ds hav e been shown to improve the scaling 
of the Nu-Ra bounds ( Constantin &: Doerind ( 1996f ): Kerswell ( 2001 )). but we shall not 
pursue such assumptions here. 

Observe that according to the above description, the trivial fields v = 0, = 
are "allowed", though they do not satisfy (|2.36p - (|2.38p for r" ^ 0. Since furthermore, 
Q^^^JAv, A6'] A2Q;^^Jv,6'] for A e M, it follows that if Q'^j^^ is bounded below, then 
the minimum is zero. 



Admissible backgrounds: 

For each i?e > (that is, for each i? > and b > 1), we thus denote a background 
field t{z) admissible if it satisfies the appropriate BCs, the same as those for T at the 
upper and lower interfaces; and if the resultant quadratic form Q'^ ^ is non-negative, 
Q;_^Jv, 9]>0 for aU allowed fields v and 9. 

Note that the positivity condition on the quadratic form Q'^ ^ [v, 6*] > is equivalent 

to 

/ Q' R [v, 611 \ 

A(r,i?o)= inf — ^p- ^ >0, (2.52) 

^-'^ , VIMI +11^11 / " 

where the infimum is taken over allowed fields v and 9. Via the associated Euler-Lagrange 
equations, the condition (j2.52p is equivalent to requiring that the linear operator 



r _ f-^S7^v + T'9 + Vp 



(2.53) 



acting on allowed fields v, 9 has a positive semi-definite spectrum, or that the lowest 
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eigenvalue Aq of Ct,r^ is non-negative. Consequently, the admissibility criterion on back- 
ground fields t{z) (for a given R c — bR/jb— 1) > 0), equivalen t to (I2.52p . is also referred 
to as a spectral constraint on t ( Doering fc ConstantinI ( 19961 )). 



Fourier formulation of admissibility condition: 

Due to the horizontal periodicity of the problem, we may reformulate the admissibility 
condition Q'^ ^ > in horizontally Fourier-transformed variables: we write the vertical 
component of velocity w and temperature fluctuation 9 as 

.Kk.'-^+Kv)yj^{^z), e{x,y,z) =y"e''-''-''+''yy^Okiz), (2.54) 



w{x,y,z) 



k 



k 



where the horizontal wave vector is k = {kx,ky) — {2T:nx/Lx,2TTny/Ly), and we write 
k'^ = |kp; we shall also write for the complex conjugate of ?ik and D = d/dz. 
We can use incompressibility to express the transformed horizontal components of ve- 
locity in terms of the vertical component, so that the admissibility criterion may be 
written completely in terms of the Fourier modes and ^k- This considerably sim- 
plifies the formulation, particularly since Q'^ is a quadratic form (equivalently, the 
Euler-Lagrange equations for the minimization problem are linear) , so that different hor- 
izontal Fourier modes decouple. The no-slip boundary condition and incompressibility 
imply that the BCs for Wii{z) are u)k — Dwk — for z = 0, 1; the BCs on 6'k obvi- 
ously depend on the choice of thermal BCs, which have so far been left unspecified. We 
note also that wq = 0; this follows from incompressibility and horizontal periodicity via 
Adw/dz = JJj^Wz dx dy = — J fj^iu^+Vy) dxdy = 0, which implies using W\z=o — that 
W = for all z. 

Substituting (|2.54p into (|2.45p and using incompressibility, as in lOtero et al. I (|2002h we 
can write the quadratic form Qt,r^ evaluated on allowed (stationary) fields v and 9 as 



^|Vvp + 2r'u;6i + |V6'|2 
Ro 



where (see IConstantin fc Doeri"^ (|l996l ); iKers"^ (|200ll) ) 

Qk = Qk;T,i?,e[f«k,6'k] = 



> 



(2.55) 



U'lu-kP + 2|Di&k|' + ;^|D'^^k|' ) + 2T'Re[wk0a 



(fc2|0kP + |D^kP 



dz ; 



(2.56) 



note that (|2.55[) is an equality for two-dimensional flows. Since the boundary terms in 

1 



(|2.49p are expressed in Fourier space as = J^k Rc[^k(-z) D^^(2:)] 

define 

to give 



z=0 



z=0 



QUJv,0] 



2t'w9+ \V9\' 



A 



> 



we thus 



(2.57) 



(2.58) 



Since the class of allowed fields [v, 9] includes fields containing a single horizontal 
Fourier mode, it is now clear that Q'^ is a positive quadratic form, Q'^ [v, 9] > for 
all allowed fields v and 9, if and only if all the quadratic forms QJ^. — Q[^.^ are positive. 
Thus the admissibility criterion for background fields t(z) (for given R^ > 0) may be 
formulated in Fourier space, as the condition that Qk[iDk,^k] > for all k and for all 
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sufficiently smooth (complex- valued) functions 'Wk(z), 9k{z) satisfying Wk — Dw^ = at 
z = 0, 1 and the appropriate boundary conditions on 9^ at z — 0, 1. 

Bounding principle: 

The thermal BCs at the plates imply an equation relating (3 and AT, whenever they 
do not specify either AT (fixed temperature) or (3 (fixed flux) directly. Thus, for instance 
(for BCs other than fixed flux), we can substitute for AT and write (|2.48p in terms of only 
/3, as done below (the fixed temperature and fixed flux cases are given in (|2.50p - (|2.5ip '). 

In such a case for a given R and b, for any admissible background field t{z) the 
inequality ^ > implies an upper bound B[t; b] on the averaged boundary flux /3. 
Via the relation between (3 and AT, this also gives a lower bound I?[r; b] on the averaged 
temperature drop across the fluid; using (|2.23p . in this way admissible backgrounds lead 
to upper bounds A/'[r; b] for Nu. For a given R> 0, the best upper bound A/'(i?) on Nu 
obtainable using this approach is now obtained by minimizing A^[r; b] over admissible 
r(z) and 50 Finally, the relationship Ra — R AT > i?P[T; b] lets us bound R, and hence 
Nu, from above as a function of Ra. 



3. Mixed (Robin) thermal boundary conditions 

The formulation and derivations above were developed independent of thermal bound- 
ary conditions, except that we have restricted our attention to fluids with thermally 
identical upper and lower boundaries, which permits the symmetry assumption t'(0) = 
t'(1) = —7. We now specialize to particular BCs to make further progress: General lin- 
ear conditions at the boundary of a fluid as in figure [T] are of mixed (Robin) type. In a 
subsequent paper we consider the more realistic case of a fluid in thermal contact with 
bounding plates of flnite thickness and conductivity. 

3.1. Fixed Biot number boundary conditions and nondimensionalization 
In dimensional terms, we choose the mixed (Robin) BCs on the plates to take the form 

T* +7fn-V*T* = a; on z* = 0, T* + ?f n • V*T* = A* on z* ^ h. (3.1) 

for some given constant < 77* < ooQ These conditions may be interpreted as Newton's 
Law of Cooling (Heating), in which the boundary heat flux is assumed proportional to 
the temperature change across the boundary: — A/n • V*T* = A/(T* — v4;*)/?7*. 

We use n = — ez, +ez on z* = O.h respectively, and nondimensionalize by substituting 
z* = hz, T* = Trcf + OT. Defining the Biot number rj — rj* /h, we fincj^ 

T-r^T,^^^^ on z = Q, T + r,T, = 'J'"' on z^l. (3.2) 

At the moment Trcf and O are still unspecified. A convenient choice, consistent with the 
nondimensionalizations introduced previously for the limiting fixed temperature and fixed 
flux cases, is to require the conducting state (u* = 0, V*T* = Ccz for some constant 
C < 0) to take the form u = 0, T=l — zin the nondimensional variables. The condition 
that T = 1 — z satisfies the nondimensional BCs (|3.2p implies {A1 — Trci)/Q = 1 + r], 

f Recall that the class of admissible r(z) depends on R through R^. 

X The limit ?;* ^ 00 is treated by writing (|3.H) in the equivalent form (for 77* > 0) 
n . v*T* -I- T* /r]* = B*^ on z* = 0, h, where (for 0< rj* < 00) = Al^/rj* . 

If There appears to be little consensus in the literature as to whether the term "Biot number" 
refers to 77 as defined in (|3.2p . or to its inverse 77"^. 
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(A* — rrof)/6 — —ri, so that (for 77 < 00) our chosen temperature scale and reference 
temperature are 

A* - a: ai + 7^{a* + a:) 

° = TW = iT^ • (^-^^ 

Having finaUy fixed a choice of dimensionless variables, the nondimensional mixed ther- 
mal boundary conditions (fixed Biot number) are 

T — ?7Tj = l + ?7 on z ~ 0, T + 7]Tz — —7] on z — \. (3-4) 

Note that the mixed (Robin) BCs (|3.4p reduce to the fixed temperature (Dirichlet) 
BCs (pTTj) in the hmit r/ 0, and to the fixed flux (Neumann) BCs (P7T^ in the Hmit 
rj ^ oo] thus we denote ry = and rj — 00 as the "fixed temperature" and "fixed flux" 
cases, respectively. 

3.2. Governing identities for fixed finite Biot number 

In the general case, neither the boundary temperature drop AT nor the flux (3 that com- 
bine in the computation (j2.23p of the Nusselt number is known a priori. However, we can 
derive a relation between them: taking horizontal averages of p.4p . and subtracting the 
upper boundary condition from the lower, we find r|z=o — T\z=i — r]lT z\z=a + T z\z=i) = 
l-|-2r;. Taking time averages and using (|2.13p and (|2.19p . we find the fundamental relation 

AT + 2r](3 = 1 + 277 (3.5) 

(this formula also holds in the fixed temperature and flux limits 77 — )■ and 77 — > cxj). 
Hence for < 77 < 00, an upper bound on (3 constitutes a lower bound on AT, and vice 
versa, and we only need to bound one of these quantities to obtain an upper bound on 
Nu = PI AT = (3/[l + 2t]{1 - P)] = AT-i + (AT"! - l) /2r]. 

Using the identity to solve for either AT = 1 + 277(1 - /3) or /3 = 1 + (1 - AT)/2r] 
and substituting into the results of Sections l2.3l - |2.5i we obtain the forms of the governing 
energy identities for mixed thermal BCs. We shall state these identities in a way that 
permits us to obtain an upper bound on /? (the relations are stated in terms of /3, 7, 
{T^\z=a,i) ^i^d (0z|z=o,i)i valid for 77 < 00)] this is the formulation suitable for small 77, 
which reduces to the corresponding previously stated identities for fixed temperature BCs 
in the limit 77 —s- 0. For completeness, in Appendix IB] we give the forms of the identities 
(equivalent for < 77 < 00) which yield the fixed flux limit 77 00. 

First, solving for AT from ()3.5p and substituting into the global kinetic energy identity 
we find 

-1^(11 Vu||2) = /? - AT = (1 + 277) (/? - 1) (3.6) 

which reduces to (j2.28p in the limit 77 ^ 0. Similarly, we can evaluate the boundary term 
in (|2.27p by using the BCs (|3.4p to solve for T at z 0, 1 in terms of T^; substituting 
into TTzl^^Q and taking horizontal and time averages, we find that the global thermal 
energy identity (|2.27p becomes 

^(11 VT||^) = {TTz\l^,) = (1 + 277)/? - v{t!U=o + t!U=i) (3.7) 

which again reduces to the appropriate fixed temperature limit (j2.29p . 

In the background flow formulation, the requirement for the background held to satisfy 
the given BCs in this case means that r(z) should obey p.4p . implying that 



At -I- 2777 = l-t- 2r; 



(3.8) 
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using ()2.32|) . Consequently, the perturbation 9 satisfies the homogeneous Robin BCs 
9 + r]n ■ V0 = at the interfaces, which in our geometry become 

9-r]9^^0 at z^O, 9 + 9^ ^ at z ^ 1, (3.9) 

and translate for horizontal Fourier modes defined in (|2.54p to the BCs 

^k(O) = 770^(0), 4(1) = -7?d4(1). (3.10) 

The boundary term in (piiS)) is thus J^9 n ■ ^9 dS = -rj /^-(n • Ve')^ dS < 0, or 

M:\1^^^-v(9!U=o + 9!U=i), (3.11) 
which can equivalently be written in Fourier space (see (|2.57p ) as 

^iLo = E Rc[4(z)D^;:(z)][^^ ^ -r^Yl (|d4(0)I' + |d4(1)P) • (3.12) 

k k 

It follows that Q't- [v, 0] > Qr.flel'^i so that a lower bound on Qt,r^ implies a lower 
bound on Q'^ the additional boundary term which appears for 77 ^ 0, cxd is stabilizing. 

Finally, the form of the governing identity (|2.48p for mixed thermal BCs, fundamental 
to the formulation of a bounding principle, may now be derived: Solving and substituting 
for AT and At using ((331) and ((3?8l) . we find 



/3AT-7AT + 7Ar = (l + 2r7)/3- 27^72. (3.13) 
Now substituting l\3.G\i and (|3.13p . (|2.48p becomes in terms of /? (compare (|2.50l) ) 

(1 + 277) (/3 -l) = b(^J^' r'2 dz-{l + 277) + 27/7^^ - ^Q^Jv, 9], (3.14) 

where we evaluate the boundary term in Q'^ _r^, [v, 9] using p. lip . 

3.3. A bounding principle for mixed thermal boundary conditions 

Following the general approach outlined in Section[2?6l the governing equations derived in 
Section [321 now allow us to formulate an upper bounding principle for the Nusselt number 
Nu in terms of the control parameter R (and hence in terms of Ra, via Ra = RAT): 

For each i? > 0, if we can choose 6 > 1 and a corresponding admissible background 
field t{z) (so that Q'^-ji^ > 0), then from p.l4p the averaged boundary temperature 
gradient /? is bounded above by 

r-l 



P<l-b+ ^ r'2 dz + 27772 j = S„[r; 6], (3.15) 

while from (|B Sp . the averaged temperature drop across the fluid AT is bounded below 

by 

AT>l + 6(2Ar-l)-6^^ f / t'^ dz + —AtA ^ VJT;b], (3.16) 

1 + 2?7 \Jo 277 J 

where the above equations define the functionals Bjj[T;b] and Vnlr^b]. (Of course, for 

< 77 < 00 the bounds (jS.lSp and (j3.16p are not independent; in principle we only 

need to find, say, the upper bound 2?,,[t;6] for /3 (for 77 < 00), since by (|3.5p we have 

Vt^[t; b] + 2r]Bn[T; b] — 1 + 2r].) An upper bound on (3 and a corresponding lower bound 

on AT then imply via (|2.23p that the Nusselt number is bounded above by 

Nu<Mr,[T;b]^Br,[T;b]/Vr,[T;b]. (3.17) 
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Slope — -7 



1(1 -Ar) 



Figure 2. The piecewise linear background profile ts{z), with r' = —7 in the boundary layer, 

and r' = in the bulk. 



Observe that for the conduction solution t(z) = 1 — z, we have ;S^[t; b] = "Drjlr; b] = 1, 
so that whenever this is an admissible profile, the bound on the Nusselt number takes 
its minimum value of 1, as expected. 



4. Piecewise linear background and elementary estimates 

As discussed in Section \2M the best upper bound J\f{R) that may be achieved by the 
above formulation for each value of R is obtained by optimizing the upper bounds on 
the Nusselt number Af[T;b] over all admissible t{z) and over & > 1. Careful numerical 
studies obtaining such optimal solutions of analogous bounding problems have been per- 
for med for plane Couette flow (which may be related to fixed temperature convection) 
by Plasting fc Kerswell ( 20031 ) and for infinite Prandtl number convection bv llerlev et al 



However, by restricting the class of admissible backgrounds t{z) over which the 



optimization is performed, upper bounds may be obtained much more readily, at the 
(likely) cost of weakening the upper bound. 

4.1. Piecewise linear background profiles 

Following Doering fc Constantin ( 19961 ) and subsequent works, we thus introduce a one- 
parameter family of piecewise linear background profiles ts{z), for which = — 7 for 
< z < S and 1 - S < z < 1: 

( Ta-j{z-S), 0<Z<S, 

t{z) = ts{z) = I Ta, 6<z<l-S, (4.1) 

[ Ta-i{z~{l-6)), l~S<z<l; 

see figure [5] Here the single parameter S < 1/2 may be interpreted as modelling the 
thickness of the thermal boundary layer. The intuition behind this definition is that 
in order for t{z) to be admissible, the indefinite term / 2t'w9 in Q'^. ji [v, 6*] (see ()2.45p . 
(|2.49p ) should be controlled by the other, positive terms. With this choice of background, 
2t'w9 vanishes in the bulk of the domain, and is nonzero only near the fluid boundaries, 
where w is small. Furthermore, since r' is piecewise constant, explicit analytical bounds 
are readily attainable, giving (non-optimal) rigorous bounds on the Nusselt number. 
From the definition (|4.ip . we immediately compute t(0) = Ta + j6, r(l) = Ta — 7(5, 
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and so 

At = t(0) - r(l) = 25-^, [ r'^ dz = 26-^^ = jAt, (4.2) 

Jo 

where it remains to choose the average = ^{t{0) + t(1)) and boundary slope 7 = 
— t'(0) = — t'(1) of the background as functions of 6. For the mixed (Robin) thermal 
BCs introduced in Section O substituting At — 25^ into the relation Ar + 2777 — 1 + 2t] 
(|3.8p we obtain the value of 7 (for given S and 77), as well as the corresponding Ar, for 
which the piecewise linear profile (|4.ip satisfies the BCs: 

1 + 2?? . (5(1 + 277) ^. 

7= \ , At = 257 = ^—-^, 4.3 

2(d + 77) + 77 

which implies the related identity 1 — Ar = 277(7 ^ 1); observe that since S < 1/2, we 
have 1 < 7 < 1/2(5 and Ar < 1. Using this 7 in the EC ([23) t(0) = 1 + 77 + 77r'(0) for 
77 < 00, we find that r^ = 1/2 (so that we may write r(0) = |(l + Ar), r(l) = i(l — Ar)), 
completing the specification of the background ts{z). 

Substituting formulas (|4.2p - (|4.3p into (|3.15p - (|3.16p and simplifying, we now find that 
the conservative bounds on (3 and AT for fixed Biot number convection with a piecewise 
linear (pwl) background profile ts take the concise forms 

/3 < Spwi,^((5, b) = B^[Ts;b] = l + b\ i-^ = i + b{j-l), (4.4) 

1 2(5 

AT > I?pwi.^((5, b) = V^[Ts;b] = l-br^ —— = 1 + b{AT - 1), (4.5) 

+ ri 

and the corresponding upper bound on the Nusselt number is Nu < A/'pwi,?)((5, &) = 
■^r][Ts]b] — Kpwi,,,(5, &)/2?pwi,r,((5, 6). Sincc 6 > 0, these bounds satisfy SpwL,,((5, &) > 1, 
1^pwi,rt{S,b) < 1, and hence Afpwi,7j{S,b) > 1, as one might expect. Observe that the 
bounds Spwi,j7(i5, b) and I?pwi,?;((5, b) do not depend explicitly on the control parameter R, 
but rather indirectly through the admissibility condition on 6. 

In the special case of fixed temperature BCs, for which AT = Ar = 1, we have 
7 = 1/2S = r'^ dz, and the bound on (3, and hence on the Nusselt number, becomes 

Nu^/3< Spwi,o(^, ^) = 1 + ^ - 1) • (4.6) 

At the opposite extreme, the fixed flux BCs impose /3 = 7 = 1, so that we must choose 
Ar — 2S — Jp r'^ dz; then the lower bound on AT (corresponding to an upper bound on 
Nu) is 

Nu-^ = AT > I?pwi,oo((5, 6) = 1 + 6(2,5 - 1). (4.7) 

Again, in this formulation the fixed temperature and fixed flux cases are the r/ — > and 
77 — > 00 limits of the bounds for general Biot number. We note, however, that the thermal 
BCs of the form (|3.4p do not specify the value of r^ in the fixed flux case rj — 00 (the 
governing equations depend only on temperature gradients, not on their absolute values); 
here we choose r^ = 1/2 for convenience, but use a different choice in the boundedness 
proof of Appendix [Xj 

4.2. Cauchy-Schwarz estimates on the quadratic form 

Recall the admissibility criterion for the background field ts{z): Q'^^ ^ [v, > for all 

allowed v and 9, or in Fourier space (by ([235l) " (|238l) ) Q'^ = Q'^,r,.R^ ["'k, ^k] > for aU 
k. For piecewise linear background fields ts(z) of the form (|4.ip . this criterion reduces 
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to a requirement that 5 is sufficiently small, for given Rc ~ bR/(b — 1). Elementary 
Cauchy-Schwarz and Young inequalities applied to the Fourier space quadratic form 
allow us to derive explicit sufficient conditions on 6 so that Q'^ > for all k, and hence 
to estimate upper bounds on Nu. 

We recall first that the boundary terms in Q'^^ ^ [v, 9] for rj E (0, oo) are nonnegative 
for Robin BCs (and vanish if = or = oo), so that it is enough to verify the 
admissibility criterion for Qr^.iie ["^; ^] (s'^e the discussion below p.lOp V Equivalently in 
Fourier space, using p.l2p to evaluate the boundary term in (|2.57p for 77 e (0,oo), we 
have 



Ql = Qi 



r;(|D0k(O)P + |D4(l)p) > Qk, 



(4.8) 



so that it suffices to obtain conditions on S to ensure Qk > for all k. To do so, we 
need to control the only indefinite term in Qk, Jq 2ri Re[wkfeli by the o ther terms. For 
completeness we review the necessary estimates from lOtero et al\ ( 20021 ): Since and 
Dt&k (and hence also Wk ^k) vanish at both boundaries, we have 



D (yik^k) dC 



< 



\wuBei\dc 



KBwul dC, 



(4.9) 



where for < z < ^ , by the Fundamental Theorem of Calculus and the Cauchy-Schwarz 
inequality we find that 



\Wkiz} 



\Bw^{z)\ 



J B^w^dC <V^^i^J \D^Wk{0\^ dCj <\/i||D2u,k 



< Vz 



|D«)k(c)r dc 



1/2 



< V^||Dti)k||j„ ij 



[0. 



(4.10) 
(4.11) 



Substituting these estimates into (|4.9p and again applying the Cauchy-Schwarz in- 
equality, for < z < ^ we obtain 



\M^)0*k{z)\ < 



< 



CdC 



1/2 



2V2 



«iiiD..kii;ij 



iD^kpdC 



?k|'dC 



1/2 



1/2- 



—\m 



a2 [0,2! 

(4.12) 



where we have also applied Young's inequality pq < ■^{ojp^ + q^/oj) for any Oj > 0. 
Proceeding similarly, we obtain an analogous estimate for i < z < 1. For the piecewise 
linear background ts{z), for which r' = —7 < for < z < (5 and 1 — S < 1, and t' — 
otherwise, applying these estimates we have 



Ts w^Ol dz 




Wk^kl dz -t- 



\w\^9l\ dz 



ai||Duik||rr 



[oa] 



1-5 

^1 



2 

[0,1] 



02 I 

/c2l 



D^wkllfoa] 



0-2 



^k||fo,l] 
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and thus 

f-i 



> 



2V2 



ai||Dwk| 



u)k^k ) dz 
1 



0-2 , 



fr2 

-I 

0-2 



(4.13) 



where norms are taken over the entire interval [0, 1] unless otherwise indicated. Substi- 
tuting this estimate on the indefinite term into Qk given by (j2.56p . we find 



2 



7(5^ ai 

2V2 

7^2 



2V2 



02 



1 - 



7(?2 02 

2\/2 

7^2 



^iiD'^kir+^fc2||^kr 



2\/2ai 



In the absence ol any additional a priori informati on, fo r instance on the deca y rate 
of the Fourier coefficients (compare iConstantin fc Doerinel ( 1996h : Kerswelll ( 200l[ )). our 
remaining estimates are necessarily fc-independent; we ensure the positivity of Qk by 
requiring all coefficients to be nonnegative. We choose ai ~ a2 ~ Il\f2\ then, using 
and dropping manifestly nonnegative terms, 



> Sk > ( ^ - 



7^ 
8 



|D^k||2+(-- 



72<5^ 



(4.14) 



We can thus guarantee that Qk > (and hence QJ^ > 0) if we choose ^'^b^ 1% < l/Rc- 
For given thermal BCs, 7 — ■y{S) is specified as a function of i5; so this is a constraint on 
6 to have Qk > 0, that is, for ts{z) to be an admissible background. Defining Sc by 



2e4 



8_ 



= 8- 



. b- 1 
bR ' 



(4.15) 



we obtain the best bound in this approach by choosing S — 6c', the piecewise linear profile 
Ts is admissible for any S < Sc- 

We observe that the estimates (|4.9p - (|4.14p . and hence the sufficient condition (|4.15p 
on S, are independent of the choice of boundary conditions on 6^ at z = 0,1 (apart from 
symmetry). However, the thermal BCs enter the admissibility condition on rg through 
the value (|4.3p of 7 = 7((5) as a function of 77. 



5. Explicit asymptotic bounds for general thermal boundary 
conditions 

Using the piecewise linear background profile ts{z) and estimates introduced in Sec- 
tion[4l we may now derive explicit analytical bounds on the growth of the Nusselt number 
Nu with the control parameter R, and hence with the Rayleigh number Ra, for thermal 
boundary conditions with fixed Biot number < 77 < 00. We begin as usual by recalling 
the results for Dirichlet (77 — 0) and Neumann (77 = 00) BCs, as in the general case it 
then becomes apparent that the fixed temperature case is a singular limit, while for any 
77 > 0, the i? — > 00 asymptotic scaling is as in the fixed flux case. In this Section we sum- 
marize the main asymptotic bounds; more detai ls including discussion of d ifferent sc aling 
regim es for < 77 < 00 will be given elsewhere (jWittenberg fc (12008!); see also I Gad 



( 20061 ) )■ while rigorous, though somewhat weaker, bounds are proved in Appendix ICl 
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5.1. Fixed temperature boundary conditions 

implies 7 = 1/2(5. 



In the case of Dirichlet BCs, we have AT = Ar = 1, i? = Ra, and 
Thus the sufficient condition (|4.15p on 5 simplifies to 5 < 5c where 



32 

Rc 



32- 



1 



bR 



(5.1) 



One c an show that the optimal choice of 5 in this formulation is 60 = 3/2 (see lWittenberg fc Gao 
(|2008[ )). for which R^ = 3R, and hence S < Sc = 4^2/3 i?-^/^ is sufhcient to obtain a 
rigorous bound. Since for this b = bo, (|4.6p becomes 



Nu^l3< 6pwi,o(<5,foo) = 1 - ^0 + ^ = -2 + 45' 



(5.2) 



for any S < Sc, the best rigorous analytical bound on the Nusselt number using this 
approach is 



13 

Nu < AApwi(i?) = Spwi,o(4, bo) = - 2 + 4^ 



_3_ 

16 



2 2 



3^6 
"32" 



Ra 



1/2 



(5.3) 

where we used the fact that for fixed temperature BCs, the control parameter R is the 
usual Rayleigh number Ra. 

5.2. Fixed flux boundary conditions 

In the opposite extreme, for Neumann BCs, we have /3 = 7 = 1, and we bound AT from 
below using ()4.7|) . Since b > 1, in order for the lower bound I?pwi,oo('5, b) = 1 — b + 25 b on 
AT to rema i n posi tive as i? ^ cx) and hence 6 ^ 0, we need b—l = 0{6). Thus following 
Otero et all (|2002r ) we choose b = 1 + cS and let c take its optimal value Coo — 1/2, so 



that (|4.7p becomes 



Nu 



AT > I?pwi,oo('5, 1 + c^S) = 1 + (1 + 5/2){2S - 1) = -S + 6^ 



1^- 



(5.4) 



The condition on S is as usual 6 < Sc, where with 7 = 1 and 6=1 + (5/2, the equation 
(|4.15p satisfied by Sc takes the form 



8 

Rc 



-R-' 



l + S/2" R 

for large R, for which (5 — > 0; and hence Sc ^ A^^^R'^^^ . Thus we have (using (|2.15p ) 

Nu-' = AT > I?pwl,oo(5c, 1 + '5,/2) ^ ^Sc ~ 

Ra = RAT > i?2?pwi,oo((5c, 1 + 4/2) - i^R^^^ 



(5.5) 



and so 



Nu < A/pwi(i?) < ^R'^' < {^R'^""^ 



(5.6) 



as m 



Otero et 



( 2002 ). Note the scaling Nu < CiR'/^ in terms of the control parameter 
R, which translates to the usual scaling Nu < C2Ra ' . 

5.3. Mixed thermal boundary conditions with fixed Biot number 

For general mixed (Robin) thermal BCs with fixed Biot number, we need to estimate 
both AT and (3, using (|4.4p and (|4.5p . where 7 and At are given in terms of 77 and S 



22 



Ralf Wittenberg 



by ()4.3|) . The sufRcient condition S < Sc for ts to be admissible, derived via the Cauchy- 
Schwarz estimates of Section IT2l is that 6c satisfies ()4.15p . which (substituting for 7 from 
(|4.3p ) here takes the form 

We shall see that in this general case with < 77 < 00, depending on the relative sizes of 
6 and 77, the scaling of the bounds behaves either as in the fixed temperature limit (for 
5 > 77) or the fixed flux limit (for 6 < rj); but that for any 77 > 0, the asymptotic scaling 
as i? — > 00 is as for fixed flux boundary conditions: 

The fixed temperature problem rj = as a singular limit: 

Recall that for Dirichlet thermal boundary conditions 77 = 0, we have AT = At = 1, 
so that we obtain an upper bound on Nu for any b > (there is no concern that the 
lower bound Ppwi.o on AT may become negative), and we can choose 6 — 1 = 0(1) for 
all 6. In this case 7; = 0, though, 7 — 1/26 is not bounded above as i? — > 00 (5 — * 0), 
and hence neither is /3; the growth in the (upper bound for) the Nusselt number in the 
fixed temperature case with increasing control parameter R = Ra is due to that of the 
(non-dimensional) boundary heat fiux. 

The situation is quite different for any nonzero Biot number 77: since < (5 < 1/2, we 
have < (1 — 26)/2{6 + 7;) = 7 — 1 < l/27y, so that now 7 is bounded above as 5 — > for 
77 > 0. On the other hand, Ar = 26^ is not bounded away from zero, so that since & > 1, 
to get a positive value for the lower bound 2?pwi,r; = 1 — & + 6Ar for sufficiently large R 
(small 5), we need b — 1 = 0{6) for each fixed rj > 0. Furthermore, we have AT 
as i? ^ 00, so that (for sufficiently large R) the growth in the Nusselt number bound 
is due to the decrease in AT, the (non-dimensional) averaged temperature drop across 
the fluid, rather than due to the growth in /3. That is, for any 7; > the (asymptotic) 
behaviour and scaling is as in the flxed flux case; the fixed temperature problem is a 
si ngular limit. (A simi lar observation was made in the context of horizontal convection 
bv lSiggers etaj] (|2004l) .) 



Scaling regimes: 

More precisely, the nature of the Nu-R scaling depends on whether (5 > 7/ or (5 < 77, 
and hence on the value of 7/: 

For sufficiently large Biot number (largely insulating boundary) 77 > 1/2, we always 
have 6 < 1]. Since for such 77, 7 is approximately constant (1 < 7 < 1 + 1/2?/ < 2; compare 
7=1 for 77 = 00), we see from (|5.7|) that a sufficient admissibility condition for Tg is 

6 < 6c = 0{Re ^^^), as in the fixed fiux case. We choose b = l+c6 < 3/2 for some c < 1, 
so P < 1 + b{'~f— 1) < 5/2 for all 6 < 1/2, and there is no transition in scaling regimes; as 
in the fixed fiux case, for all sufficiently large 77 the growth in Nu is due to the decrease 
in AT. 

For relatively small Biot number (largely conducting boundary) 77 < 1/2, on the other 
hand, it is possible to have 6 > rj for low enough thermal driving, and thus distinct 
regimes exist. In particular, consider the case of small Biot number (77 <C 1, near the 
fixed temperature limit), where we can identify two distinct scaling behaviours: 

• "Fixed temperature scaling": For sufficiently small R, we have 6 ^ rj, so that 7 ^ 
1/25 and At '--^ l|3 and the sufficiency condition (|5.7p is 6 < 6c ^ 0{Rc^^^). Since At is 

f Proceeding more carefully, for 5 > 77, we have 1/45 < 7 = (1 + 2t])/2{5 + 7;) < 1/5 and 
1/2 < Ar = (5(1 + 2r))/{5 -h 77) < 1 
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bounded below away from zero, so is the lower bound I?pwi.?)('5, b) — H-6(At— 1) > 1 — 6/2 
on AT for any fixed b < 2. Thus we may obtain a bound on Nu in this regime by choosing 
any & g (1; 2), and by comparison with the fixed temperature problem, it is sufficient to 
choose b-l = 0(1), in which case we have Sc = 0{[{b - l)/bR]^/^) = 0{R-^/^). While 
AT = 0(1) (so that Ra = 0(i?)), we have that /3 < i3pwi,,)(<5, 6) = 1 + 6(7 - 1) = 0{b/d) 
grows as (5^^. Thus clearly when 77 <C 1, for sufficiently small but increasing R, the scaling 
properties are as in the fixed temperature case, and the growth in Nu = /3/AT is driven 
by that of (3. 

As the driving R increases, 5 decreases, and eventually becomes less than the Biot number 
rj; based on the fixed temperature scaling 5 = 0{R~^^'^) = 0{Ra~^^'^) the transition at 
S = r] occurs when 77 = 0{Ra~^^^), or Ra = 0{r]~^). 

• "Fixed flux scaling" : Once the "boundary layer thickness" 5 has decreased below 77 > 
for increasing i?, we enter another regime (which does not exist in the fixed temperature 
case 77 = 0), in which for fixed rj the growth in 7 saturates, while Ar — 0{5) decreases. 
Asymptotically for (5 ^ 77, we have 7 ~ (1 + 2ri)/2rj = 7max('7)j while Ar ^ (5(1 + 2?7)/77, 
and for each fixed 7; > the behaviour is now as if we had Neumann thermal BCsyJ 
More generally, for < 77 < 1/2 and decreasing 5 < 77, we have 7 = 0{ri~^) and 
Ar — 0{5/if). In order for the lower bound Ppwi,?) = 1 — 6+5Ar on AT to remain positive 
as (5 ^ 0, we must choose 6=1 + 0{S/ri), so that /3 < Spwi,,, = 1 + 6(7 — 1) = 0(77"^) 
saturates, while AT > 0{S/ri); hence the growth in Nu is now due to the decay in 
AT, as in the fixed flux case. In this regime the scaling behaviours are Ra > 0{SR/ri), 
i?e = 0{rjR/S) and S = 0{j-^/^ R7^^'^) = 0{r]^/^R-^/^) = 0{Ra-^/^); more precise 
asymptotic statements are given below, with weaker, but rigorous results in Appendix [Cl 



Asymptotic scaling of bounds for < 7/ < 00; 

Having outlined the behaviour in the different regimes, we here derive the scaling of 
the bound on the Nusselt number in the limit of large driving, R 00, so that (5 <C 1 
and i5 <JC 77, deferring a more detailed discussion of scaling behaviour in the different 
regimes using this Cauchy-Schwarz analysis, and a com parison with numerical solutions 
for piecewise linear backgrounds, to IWittenberg fc Gao (^008). 

In the light of the above discussion, for (5 ^ 7y we must choose 6 = 1 + cS, where the 
optimal value of c turns out to be 

27? 



1 



" 477 ■ 

Using this optimal choice of 6, the lower bound (|4.5p on AT becomes 

<5(1 + 277) 



(5.8) 



AT > Ppwi,^((5, 1 + c„(5) = -c„(5 + (1 + CrjS)- 
_ 5(1 + 277)3 + 2,5 



5 + 77 



(5 + 77 
3 (5(1 
' 4 



277) 



(5.9) 



while similarly, the upper bound (j4.4p is 

1 + 27? 



Crjd + (1 + C^(5) 



1 + 277 

2((5 + 77) 



2((5 + 77) 



1 + 277 
277 ' 



f More precisely, for 5 < rj, we have "fma 
57max = 5(1 + 2r;)/277 < Ar < (5(1 + 277)/77 : 



,/2 = (1 + 2r?)/4r; < 7 < (1 + 2r/)/2r; = 7„ 

= 2(57inax. 



(5.10) 



and 
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so that an upper bound on the Nusseh number for admissible S t] is 



compare (|5.2p and ([57 

Observe that the width 5bl of t he thermal boundary layer is often related to the 
Nusselt number via 5bl = {2Nu)~^ ( Niemela &: Sreenivasan our high-i? result 

for the piecewise linear background, S ~ {3Nu/2)~^ for 77 > (or S ^ {4:Nu/3)~^ for 
77 = 0), may be interpreted as a systematic statement of such a boundary layer model. 

Returning to the computation of asymptotic bounds, we note from (|5.8p that for r] > 
1/2, = 1/2 + 1/477 < 1: while for rj < 1/2, c^J = {l + 2r])5/Ar] < 6/2t], so that whenever 
S <C min(77, 1) we have c^(5 ^ 1; consequently 6 = 1 + c,,,(5 ~ 1 and = h R/{h — \) ^ 
R/cjjS. In this case the condition (|5.7p is thus 

'^^ - 32#±?4i?-i ^ 32^i^i±^5i?- = 8^5i?-, (5.12) 
(l + 2r/)2 ° (1 + 277)2 477 1 + 277 

or 6c ~ 277^/'^(l + 277) ^^/"^i?^^/"^. Substituting into the above bounds, we have 

i?a = RAT > i?Ppwi,„('5e, 1 + c,,Sc) ^ ^i±^ J^i? ~ ^ ('i±^) ^ i?2/3^ (5^^4) 

4 77 2 \ T] J 

SO that we obtain a bound on the asymptotic scaling as i? ^ cxi of the Nusselt number 
with the Rayleigh number whenever rj > 0: 

Nu < I ('-±^) ,/| ^'^ Ra-I^ = JlRa-l\ (5.15) 



3 V / V 3 V 1 + 2?7/ V 27 

independent of the Biot number. Observe in particular, by comparison with (|5.6p . that 
the prcfactor ^2/27 is the same as for the fixed flux problem. 



6. Conclusions 

In formulating the energy identities and bounding problem for the Rayleigh-Benard 
model with finite Prandtl number and general thermal BCs at the upper and lower 
boundaries of the fluid, we have demonstrated that the fixed temperature and fixed flux 
extremes may indeed be treated as special cases of a more general model, within which 
one can rigorously prove energy boundedness and bounds on convective heat transport, 
and obtain asymptotic scaling results; we expect that such an approach may be applicable 
to other related convection problems. 

While the scaling of these analytical bounds on the Nu-Ra relationship remains well 
above that observed experimentally or in direct numerical simulations, some of the qual- 
itative conclusions may be instructive. Of particular interest is that — at least for the 
piecewise linear backgrounds ts{z) treated here — while for each fixed R the bounds de- 
pend smoothly on 77 for < 77 < 00, the asymptotic R —^ 00 scaling of the bound for any 
nonzero Biot number is as for the rj = 00 fixed fiux problem. That is, the limits 77 ^ 
and i? ^ 00 do not commute: fixed temperature conditions in fact form a singular limit. 

Furthermore, the bounding calculation indicates the existence of two distinct scaling 
behaviours for sufficiently small nonzero rj; it would be of interest to observe these in 
fixed Biot number direct numerical simulations: For small Rayleigh number iZa, there is 
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a "fixed temperature scaUng regime" in which the usual assumption of Dirichlet thermal 
BCs is approximately valid, as the growth in the convective heat transport measured by 
Nu is largely due to the increase in the averaged boundary heat flux /3. As the control 
parameter i?, and hence i?a, increases, a transition occurs when the "boundary layer 
width" 5 becomes comparable to rj (in our calculations this occurs for Ra — 0{ri~'^)), 
beyond which a "fixed flux scaling regime" is entered, in which further increases in 
Nu are driven by decreases in the averaged temperature drop AT. This observation 
provides mathematical support for the heuristic argument that when the Nusselt number 
is sufficiently high, the boundaries act effe ctively as insulators. 

In the three-dimensional simulations of ' Verzicco fc Sreenivasan ( 2008 ) in cylindrical 
geometry with perfectly insulating sidewalk and a perfectly conducting upper bound- 
ary, replacing the lower fixed temperature BCs with fixed flux conditions was observed 
to have little effect on the heat transport for a given Rayleigh number for sufficiently 
small Ra < 10^, and to decrease the transport fo r Ra > 10^. In cons t rast, the two- 
dimensional, horizontally periodic computations of Johnston fc Doering ( 2008[) showed 
essentially identical heat transport for fixed temperature and fixed fiux BCs at both up- 
per and lower plates for 10^ < Ra < lO^*'. In this context we observe that the prefactor 
in our asymptotic analytical bound Nu < C Ra ' increases from Co = 3%/6/32 w 0.230 
to Cri — Coo — \/2/27 « 0.272 for rj > 0; that is, within the framework of our upper 
bounding calculations with piecewise linear background it appears that the estimates on 
the heat transport increase when the boundaries are not perfectly conducting. It remains 
to determine whether this increase is an artifact of the choice of background t(z) or of 
the background flow bounding approach in general. A further consideration is how the 
presence of finite width conducting plates (see Part 2 of this work) modifies conclusions 
obtained with the fixed Biot number simplification. 
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Appendix A. Boundedness of ||T|p and ||up 

For completeness of the rigorous argument, we show the uniform boundedness of the 
temperature and velocity fields, which we may state as a theorem: 

Theorem 1. The velocity field u and temperature field T satisfying (|2.7p - (|2.9p for 
finite Prandtl number, < Pr < oo, and for no-slip velocity boundary conditions and 
thermal boundary conditions of general Biot number < r] < oo, are uniformly bounded 
in L^. 

Remark: Such boundedness has already been shown for Raylcigh-Benard convection 
with fixe d temperature BCs by K erswcU (2001), following th e unde rlying approach intro- 
duced by lDoering fc ConstantinI |l992h (based on an idea of lHopj (jl94lh ) in the context 
of shear flow. However, in both of these cases the Dirichlet boundary conditions allow 
ready control of the indefinite term in real space; since for general thermal BCs we are not 
assured control of 9 at the fluid boundaries, in our proof instead w e use incompressibility 
and Fourier space estimates based on those of lOtero et ai\ (|2002h . 



Proof. We begin the demonstration of Theorem [T] by reviewing the basic problem 
formulation and identities: With thermal boundary conditions imposed at the upper 
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and lower limits of the fluid, we consider horizontally periodic temperature and velocity 
fields T(x, t) and u(x, i) satisfying ()2.7p - (|2.9|) . where u satisfies incompressibility and 
no-slip BCs. Choosing a background temperature field r(z) which satisfies the given 
thermal boundary conditions, we define v(x, t) and 0{x,t) via the decomposition (|2.33p . 
u(x, t) = v(x, t), r(x, t) = t{z)+9{x, t), and thus obtain the evolution equations (|2.39p - 
(|2.40p for their norms: 

^ '|v||' = -4l|Vv||V / w9, (Af) 



2PrRdt" " R 

1 d „,„2 '\V9f +A0e:\'- 1 0,y + A0T'\l_~ I wOt'. (A2) 



2dt 1^=0 y 1^=0 

We form the linear combination (lA 1\ + /i- (jA 2\ . where the weight /i will be chosen later: 

f d 
2di 



^ II l|2 
V 



PrR 



R 



^ II Vv||^ + fi w0- t'wO - \\vef 



T'e,^ + Aee,\l^^ + AeT'\l^^. (A3) 



We choose, as before, a piecewise linear r(z), defined as in (|4.f p and figure [5] for 
< (5 < f/2: 

r T,~7(z-,5), 0<z<5, 
t(z) = < Ta, (5 < z < f - 5, (A 4) 

[ T, -7(z-f + ,5), f-(5<z< f; 

As in Section SjI we find using the thermal BCs that (143]) 7 -t'(0) = -T'(f) = 
(f + 2rj)/2{5 + ry), while for r/ < 00 we also have Tq = 1/2; for now we defer the (at 
present arbitrary) choice of Ta in the fixed flux case r/ = 00. 

A.l. Estimates independent of thermal BCs: 

The estimates on the indefinite quadratic terms are performed in Fourier space using the 
definition (|2.54p . while other terms are readily controlled in real space; thus we split the 
dissipative terms as 

l|V0ir= Q + j \V6\^ ^\\\ve\\' + \a^ j'^ (k^\e^\^ + \Dk\^) dz (A 5) 

and, using incompressibility, 

||Vv||2>i||Vvf + 1^^] (fc2|z;;k|2 + 2|Du;kP + ^|D2^«kp) dz, (A 6) 

with equality for two-dimensional flows. 

We bound / wO using estimates of the form (|4.12p . which imply that for k ^ 0, any 
< (5 < f /2 and any p, g > 0, 



/ \w)^9l\dz+ / \w)^ei\dz 

, Jo Jl-5 , 



6^ 

< 



4^/2 



p\\Bwa' + -||d4|P + i^llD^iikf + -ll^kf 
p k'^ q 



(A 7) 
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In particular, using (|A 1\ with 5 = 1/2, recalling wo — 0, 

^ J ^ Jo ^ Jo 



< 



-y 



16V2 



aillDiikir + — ||D0k 
ai 



2 , ||T^2 |i2 , ^ 



fc2 



02 



Similarly, using the definition (|A 4p of t{z), as in (I4.13P we have 



1 
A 



t'wO 



< 



IE 



< 



2 

k 

4^2 



E 



T'(l()k6'k + Wk^k) rf^ 



flsllDiikf + -IID^IP + ^WB^Wkf + -||0k| 

03 04 
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(A 8) 



(A 9) 



The other two terms in (jA 3p containing r' may be estimated directly, since for any 
05 > 0, we have (also using II^^H < ||V0||) 



T'e, + AT'0\^^^^jA{ 0,dz+ e,dz\--iA e,dz 



1-5 



-7A 



l-<5 



1-5 



< 7^1/Vl - 25 eldxdydz 

<iA^'^\\e,\\ < I (a5-f'A + -\\Ve\\' 



9z dx dy dz 
1/2 



(A 10) 



Substituting (|A 5p - (|A lOp into (|A 3p . and collecting like terms, we thus obtain 



ld_ r 
2di 



ll^llV^llvir 



< o«57'^ - Trt; II Vv|r - ( 1 - - ) II V0|r + A 



2R 
A* ,2 



05 



lz=0 



/i a2/i 047^2 \ 1 



Y 16\/2 4V2 ; P 

^Vfi-^^ ^ 

^ V2 16V2a2 4\/2 04 

7^2 



2 16\/2ai 4^/2 



-llD^^kll 

fc'ii^kir 

llD^wll^ 



03 



(All) 



At this point we are free to choose the constants 01-05, as well as fi and S. It is convenient 
to begin by selecting 01-04 so that /x/16\/2oi.2 — 7(^2 /4\/2 03,4 = 1/4, and then, after 
substitution, to choose fj, and S to satisfy a2/i/16\/2 = a4'yS'^ / 4^/2 — ii/iR. This gives 



oi = 02 = 



4^/2' "^-"^-^2' 
and 72(5''/8 = ^/4i? = 8/i?^^ so that 7(5^ = 8/R. Choosing, furthermore, 05 — 2, and 



752 
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substituting, (jA becomes 
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1 d 
2di 



32 n 1,2 
V 



<l'A^^\\V^rf-\\\V9\\' + Ae94l^^ 



iDwkl 



<7'v4--glivvf -IiIv^IIVa: 



(A 13) 



A. 2. Poincare and related inequalities 

It remains to establish Poincare- like inequalities controlling ||v||^ and H^H^ by ||Vv||^ and 
II V^ll , respectively. In the following we consider only functions periodic in the horizontal 
directions, on {x, y) S [0, L^] x [0, Ly], so that further discussion of "boundary conditions" 
refers to the vertical boundaries at z = and z = 1. 

For nonzero functions "0 which vanish at the vertical boundaries, we have the Poincare 
inequality || VV'|P/||V'lP ^ -^d, where Ad — tt^ is the lowest eigenvaluc0 of — on this 
domain [0, L^,] x [0, Lj,] x [0,1] with homogeneous Dirichlet BCs sX z — 0,1. Applying 
this inequality to the three components of the velocity field with no-slip BCs, we have 



ivvir < -A 



(A 14) 



In order to establish an analogous inequality for — ||V6'|| , and to find 7 and 6 from 
7(5^ = 8/R, we need to consider the different thermal boundary conditions separately: 



Inequalities for fixed temperature BCs: 

For functions 9 satisfying Dirichlet BCs 6*12=0,1 = we have, as discussed above. 



\V9\\' < -Xd\ 
1 1 



(A 15) 



for Xd = tt^. Since in this case the term A 9t'\___^ in (|A 3|) vanishes, we can improve 
upon (|A 10|) to give, for < (5 < 1/2, 



9z dx dy dz 



1 

1-5 



< 7AI/2JI/2 



9z dx dy dz 

1/2 



9l dx dy dz 



9l dx dy dz 



i-s J J A 



1/2 



1 



1 



< - 2a57^M+ — \\V9 



05 



(A 16) 



Using (|A 16p instead of (|A lOp . we obtain an equation similar to (|A lip in which the first 
term on the right-hand side is a^^^'^SA, and choose 01^05, /i and 5 as before. In the fixed 
temperature case we have 7 — 1/25, so that the condition on 5 becomes 7(5^ — S/2 — 8/ R, 



t The use of A (or Xd,n,r) to represent eigenvalues only in this Appendix |X] should not be 
confused with the conductivity A elsewhere. 
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01 5 — 16/R, 7 = R/32. Hence the equivalent of (lA 13|) becomes 
32 „ 



1 d 
2di 



m 



PrR^ 



1 

z=0 













32^ 






< 






32 






< 


^A 




32^ 



Xd 
4 

4 



32 



where we used 
((XT4)) ^ (IXT5l ). 



,1 

lz=0 
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(A 17) 



for Dirichlet thermal BCs, and the Poincare inequalities 



= 0. In addition, 



Inequalities for fixed flux BCs: 

When 9 satisfies Neumann BCs 6z\z=o,i = 0, we again have 
7=1, and 5 is found from •jS'^ = = 8/R, so S — hence (lA 13p becomes 

n2 32 „ 



ld_ 
2di 



PrR^ 



<A 



(A 18) 



In this case, while the Poincare inequality (|A 14[) holds for the velocity field as before, 
the temperature field requires a bit more care, since under Neumann thermal BCs the 
equations are invariant under 9 i— > 0+constant, and we have no immediate control of 
ll^ll" by ||V0|| . However, in general for nonzero functions ip with mean zero, we have 
II VV'II^/IIV'II^ ^ ^N, where Ajv = tt^ is the lowest nonzero eigenvalue of — on this 
domain with homogeneous Neumann BCs at z = 0, 1; and we can satisfy the additional 
condition on the mean (in the fixed flux case only) by exploiting the remaining freedom 
in the definition (|A 4p of the background r(z): 

Since for fixed flux thermal BCs the flux out at the top of the fluid exactly balances 
the flux in at the bottom, the total heat content is preserved over time; more precisely, 
letting Ta — A~^ J T he the mean temperature over the fluid, in this case (|2.18p becomes 
dTa/dt = 0. Thus we may choose the (previously arbitrary) average of t(z) to be the 
(constant) average temperature, A'^^ J t = Ta = Ta, thereby completing the definition 
(jA 41) for 77 = cxD. Hence by construction the perturbation 9 has mean zero, J 9 = 0, so 
that we have the inequality 



Substituting (|A 14|) and (|A 19|) into (|A 18|) now gives 



(A 19) 



ld_ 

2dt 



32 



<A 



<A 



1 

4 



XN\\ef + ^Xn\W\ 



2Pr 



i?2 
Xd 



X 



N 



32 



(A 20) 



Inequalities for mixed thermal BCs: 

For thermal BCs with general Biot number < 77 < 00, we do not have such simple 
expressions for 8 and 7, but as before we define 5 = 6{R, rf) by 7(5)(5^ = (1 + 2ri)5'^ /2{5 + 
rj) = 8/R. Weakening the inequality (|A 13^ by using an upper bound for 7 (valid for 
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77 > 0) uniform in (5 > 0, we have in this case 



1 d 
2dt 



32 
PrR^ 



^ {^ + 2v)\ ^ 
- 4,^2 



A 



(A21) 



Observe that for 77 S (0, 00), the boundary term in (jA 2ip does not vanish in general; in 
heu of a Poincare inequahty, we exploit this fact in the following to control 

For < 77 < cxD, it is straightforward to verify that for nonzero functions 7/; on a general 
domain il, the stationary values of the ratio 



A = 



are the eigenvalues of — on Vl with homogeneous Robin BCs on dVl, and are attained 
at the corresponding eigenfunctions; that is, at the solutions of 

-V^t/j = At/; in fi, n • Vt/; + 77" V = on dO., 

so that the lowest eigenvalue Xr^q, > minimizes the above ratio. Specializing to our 
fluid domain [0, L^] x [0, Ly] x [0, 1] and evaluating the boundary term, we have similarly 
that for nonzero horizontally periodic functions ip, we have 

J\V4>\^ + Tj-^ J^4>^ds _ 11^11' + 77-1^(^^1,^0 + 



> Xf 



where Aij > is the lowest eigenvalue of — on this domain with homogeneous 
Robin BCs n • Vt/j + r]~^4' = at z = 0, 1, or equivalently ip — r]^z = at z = 0, 
ip + vi'z = at z = 1 (see p.9|) ). In particular, the temperature perturbation field 
6 for our convection problem with mixed thermal BCs satisfies Aij||6'|p < ||V6'|p + 

rj^^A (j)'^\z=Q + 0^\z=i^ ■ However, since itself satisfies the Robin BCs (|3.9p . we have by 



(IB51) that < v^^a(^9^\z^o 
BCs we have the inequality 



= ~A 



\z=0 



< -X 



^; consequently for mixed thermal 
f. (A22) 



Substituting (|A 14p and (|A22p into (|A2ip . we thus estimate 



ld_ 
2dt 



32 



< 



< 



(1 + 277)^ 
4772 

(1 + 277)2 
4^2 



A- 
A- 



1 

4 



II ||2' 



min I 2Fr^,l 
4 \ Xr' 



32 



z=0 
||2 



(A 23) 



Uniform boundedness for general thermal boundary conditions applied to fluid: 

The estimates (|A 17p . (|A 20p and (jA 23P are all differential inequalities of the form 
dE/dt < 61 - b2E for constants 61,62 > 0, where E{t) = \\e\f + (32/Fri?2) It 
follows, using GronwaU's inequality, that E{t) < E{0)e~'"'* + (6i/62)(l - e"''"*), thus 
completing the proof that v and 9, and hence u = v and T — t + 9, are uniformly 
bounded (by E{0) + 61/62) in L2 for aU t. □ 
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Appendix B. Governing identities valid in fixed flux limit 

In Section 13.21 the energy identities required to formulate a bounding principle were 
stated in a form valid for < < cxd, which are appropriate for exploring the fixed tem- 
perature limit rj ^ 0, and let us study the effect of finite conductivity as a perturbation 
of the ideal case of perfectly conducting boundaries. By instead writing these governing 
identities in terms of AT, At, (r^|z=o,i) and {0'^\z=o,l)^ we obtain a formulation valid 
for Q < rj < oo and relevant to the insulating boundary fixed flux limit — > oo. For com- 
pleteness, we state these identities here; of course for < < oo the formulas (|B ip - (IB 8|) 
in terms of AT are equivalent to (|3.6p - p.l4p in terms of /3. 

Solving for (3 from (|3.5p and substituting into the global kinetic energy identity (|2.25p , 
we first obtain 

-^(11 Vu||2) = /3 _ AT = ^(1 + 2r,) (1 - AT) , (B 1) 



,1 



which reduces to (|2.30p in the limit 77 — > oo. Correspondingly, evaluating {TTz\^_^ by 
solving for Tj at z = 0, 1 in terms of T and averaging, the general form of the global 
thermal energy identity (|2.27p which reduces to the fixed flux limit (|2.3ip is 

i(||VT||2> {TTz\\j = AT+ ^ (1 + AT) - i(T^U=o + T^U=i). (B2) 

The homogeneous Robin thermal BCs satisfied by the fluctuation 9 in the background 
flow formulation (with a background field r(z) obeying the identity (|3.8p ) may be written 
for 7? > as n • W6 + 77^^^? = 0, or in our horizontally periodic geometry 

e,-r]-^9 = at z = 0, 6*^+ 77-161 = at z = 1; (B3) 

the individual horizontal thermal Fourier modes thus satisfy 

D^k(O) = 77-1 4(0), d4(1) - -rr' 4(1)- (B 4) 

The appropriate formulation of the boundary term in (I2.49P valid for < 77 < 00 is thus 
Jj9n- V0 dS = —11^^ Jj 0^ dS < 0, which in our geometry becomes 

^l!-o = --(^l-o+^l-i) (B5) 

IZ-U ^ 



^E(l^k(0)|^ + |4(l)p) <0 (B6) 



^ k 

in real and Fourier space, respectively, again verifying the stabilizing effect of the bound- 
ary term in Q'^ Substituting for P and 7, we may write the identity (|3.13p instead in 
terms of AT and At, as 

/3At - 7AT -f 7AT = —(1 + 277) (2 At - AT) - — At^. (B7) 
277 277 

For mixed (Robin) thermal BCs with constant Biot number 77 > 0, we can now substitute 
(|B ip and (IB 7P to write the governing identity p.48p in terms of AT as 

J2 J., , --r 277 , 



i-(l + 2,)(l-AT)=5(£ 



(compare (I2.5ip ). using (|B Sp to evaluate the boundary term in Q'^ _r„[v, 0], which allows 
us to derive the lower bound (|3.16p on AT; of course for < 77 < 00, this expression is 
completely equivalent to p.l4p . 
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Appendix C. Rigorous admissibility conditions and bounds 

In Section [5l explicit analytical bounds on the dependence of the Nusselt number 
Nu on the control parameter R and hence on the Rayleigh number Ra were obtained, 
using a piecewise linear background ts{z) and elementary estimates on the quadratic 
form. However, for fixed flux and general Biot number thermal BCs (77 > 0), the condi- 
tions (|5.5p , (|5.12p on 6 and the bounds (|5.6p , (|5.15p were derived using some asymptotic 
approximations for 6 <^ min(r7, 1) (that is, for sufficiently large i?), and are thus not rig- 
orously applicable for fixed nonzero 5. The arguments of Section [5] may however readily 
be adapted to give rigorously valid admissibility conditions on ts, and hence rigorous 
bounds on Nu{Ra), at the cost of weakening the prefactors. 

We choose to present the rigorous bounds as results valid uniformly in 77 and for 
< (5 < 1/2; the prefactors may be improved by 0(1) corrections in several places by 
restricting the range of rj and/or S values. In our formal development we account for all 
relevant scaling regimes by using a balance parameter b = bc{r], 5) defined on 5 S (0, 1/2], 
?7 e [0, 00] via 

{\ forO<77<i, ^>77, 

6c-l = < 1^ forO<77< i, ,5<r7, (CI) 
I ^ for -i < 77 < CX3 (so necessarily ^ < 77) . 

Note that the function bc{rj,5) is continuous on its domain and agrees with the values 
used previously in the limiting fixed temperature and fixed flux cases, 60 3/2 and 
''oo = 1 + Ccyo5 = 1 + 5/2, respectively. We immediately observe that 

l<fcc<| (C2) 

for 77 > 1/2 this follows from (1 + 277)/4?7 = 1/2 + I/477 < 1 and (5 < 1/2. 

In the first lemma we estimate the formulas for upper bounds on Nu and lower bounds 
on AT for piecewise linear backgrounds uniformly in terms of 5 and 77: 

Lemma 2. For general mixed thermal BCs with arbitrary Biot number rj g [0, 00], the 
upper bound Afpwi,ri on the Nusselt number and lower bound Ppwi,?; on the temperature 
drop across the fluid, obtained using a piecewise linear background profile of the form ts 
(|4.ip with 6 G (0, 1/2] and a balance parameter b — bcijj, 5) defined by (|C 1[) . satisfy 

AApwu,(<5,6c) < ^, 2?pwi,,(<5,6e)> ^^^^^^1^. (C3) 



Proof. For piecewise linear backgrounds T5, with 7 and Ar given by (|4.3|) . the upper 
bounds on (3 and AT from (I4.4p - ()4.5|) using b = be satisfy 



^pwi,^(<5, 6c) = 7+(^'c-l)(7-l) = 



1 



2{S + rj) 



[l + 2r;+(fee-l)(l-2<5)] 



s 

2?pwi,,(<5, 6c) = At + (6c - 1)(At - 1) = —— 1 + 277 + (6c - 1)^(25 - 1) 

+ 7J I 



SO that 



A/'pwKr,('5, 6c) 



Bpwi.r,(<5, 6c) _ 1 1 + 27/ + (6c - 1)(1 - 25) 



2?pwL^(<5,6c) 26 1 + 2?7+ (6c- 1)(25- 1)77/5- 
Now we consider the three cases in the definition of 6c: 
I: For < 7? < 5 < 1/2: 
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We have 6c - 1 = 1/2, so 1 + 2?; + (6c - 1)(1 - 2(5) = 3/2 + 2?? - 5 < 3/2 + 2-q and 
1 + 277 + (6c - 1)(2^- 1)??/^ = 1 + 3?7- ?7/25 > 1/2 + 3??, so that 



A/'pwi,,,((5, 6c) < ^ TT"^ < T^, 2?pwi,r,('5, 6c) > 



2(5 1 + 67; -25' -pwM/v",".. - 



37? 



1 5 1 
- 2 (5 + 7) - 4' 



II: For < 5 < ?7 < 1/2: 
We have 6c - 1 = 5/2ri, so 1 + 277 + (6c - 1)(1 - 2(5) < 1 + 277 + 6/2ri < 3/2 + 277 and 
1 + 277 + (6c - 1)(2(5 - 1)77/(5 1/2 + 277 + (5 > 1/2 + 277, so that 



^ + 277 



1 S S 

> - T > • 

- 2 (5 + 7/ - 477 



AAp„i,,((5, 6c) < ^ ^ < ^, I?pwi,,(5, 6c) > 

III: For < (5 < 1/2, 77 > 1/2: 
We have 6c - 1 = (1 + 2T])d/4r] = c^6 (see dl^)), so 1 + 277 + (6c - 1)(1 - 2(5) = 
(1+277) [l+(5/477(l-2(5)] < (l+2?7)[l+(5/47/] < (1+277)5/4 and l+277+(6c-l)(2(5-l)77/(5 = 
(1 + 2?7)[1 + (25- l)/4] = (1 + 277) [3/4 + (5/2] > (1 + 277)3/4, so that 

14 + (5/t7^ (53 31+ 277 

K^UW < < ^, > ^4(1 + 277) > -^S. (C4) 

□ 

Remark: The form of the bounds in (jC 3p is chosen so that A/'pwi.r; is uniform in 77, and 
I'pwi,?) is continuous on 77 e [0,oo], (5 G (0, 1/2]; however, the prefactors Ci in the bounds 
of the form A/'pwi.r; < Ci/5, 2?pwi,?) > C'2(6c — 1) may certainly be improved for particular 
values or over restricted ranges of 77 by 0(1) factors from their values Ci — 3/2, C2 = 1/2 
in (IC3l) . 

For instance, for fixed temperature BCs ?7 = 0, from (|5.2p we have Afpw\,o < 3/45 and 
'Z'pwi.o = 1; while for fixed flux BCs 77 = 00 we can show from (I5.4p that A/'pwi,oo < 2/35 
and I?pwi,oo > 35/2. Furthermore, case III in the proof of the above Lemma shows clearly 
in (|C4l) that for 7/ > 1/2, it is sufficient to take Ci = 5/6, C2 = 3/2. 

Lemma [2] gives an upper bound Afpwi,ri on Nu and a lower bound Ppwi.?) on AT in 
terms of 5, provided 5 is such that the piecewise linear background ts is admissible for 
i?o = 6i?/ (6—1) where 6 = 6c is defined in (|C ip . The next result gives sufficient conditions 
on 5, of the form 5 < ^s, for admissibility of ts as a function of R and 77, using a balance 
parameter 6s which coincides, at 5 = 5s, with 6c used in Lemma [2j 

Lemma 3. Consider Rayleigh-Benard convection subject to thermal boundary condi- 
tions with Biot number rj G [0, 00] . 

(a) For each Re > 0, the piecewise linear background field ts{z) defined in (14. ip is 
admissible, Q'^ ['^: ^] ^ for all allowed fields v and (see Section \2.6\) . provided 
5 < min(5c,l/2), where Sc satisfies (15. 7^ ; 



= ~ (1 + 27^)2^'= ■ 

In particular, for any R > and b > I we may choose R^ = bR/ (b — 1). 

(b) For each 77 > and R > 0, a sufficient condition for the piecewise linear background 
Ts to be admissible for R^. 6si?/(6s — 1) is that 5 < min(5s,l/2), where bs{rj,R) and 
ds{r/,R) are defined as follows: 

I: ForO<r]< 1/2 a7id R < Rsirj) = (8/3)7/-^ (where we define R,{0) = 00): 
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II: For 0<T]< 1/2 and R> Rs{r]) = (8/3)r/-2; 



///; For 1/2 < 77 < oo and all R: 

Furthermore, for < 1] < 1/2 we have that Ss{ri,R) is continuous in R, and 6s > rj if 
and only if R < Rsifj)- 



Proof. Part (a) of the lemma was proved in Section 14.21 using Cauehy-Schwarz esti- 
mates on the indefinite term in Qk, where we also used the relationship from Sect ion [2. 6 1 
between admissibility (Q'^ ^ > 0) and positivity of the quadratic forms in Fourier space, 
Qk > for all k. 

To prove part (b), we first observe for < < 1/2 that for R — i?s('7), we have 
2(2/3)1/2^7^^^ = = 2{r]/3y/^R^^^^; hence (5s(ry, R) is a continuous function of R, and 
Ss> ri when R < i?s(?7), while for R > Rs{r]) we have 6s < rj- The function 6s thus agrees 
with bc{r},6s) from (|C ip . and as in (|C 2p we have that 1 < 6s < 3/2. 

Next, we define Re{r], R) = bsR/ {bs — 1) as in the statement of the lemma, and then let 
6c = 6c{r], R) be the critical value of 6 as in part (a) for this Re = Reiv^ Then since the 
result of (a) indicates that Tg is admissible for this Re whenever 6 < mm{6c{ri, R), 1/2), 
to conclude the proof of part (b) it is sufficient to show that 6s < 6c in each of the three 
cases: 

I: For < 77 < 1/2 and R < Rs{ri): 

II: For < 77 < 1/2 and R > Rs{r]): 

(1 + 2r/)2 6s - 4 2ry 3 3 ' " 

III: For 1/2 < 77 < c» and i? > 0: 

. (5. +77)2 (1 + 277) 5s 1 

'(l + 2r/)2 477 6s -"(1 + 277)2 77 3" 31 + 277' 



6f. ^ 32 (i±Mi.ii?- > 8.-4^ (i±M^?i?- ^ ^-^i?- 5s = 6t 

□ 



We may now combine the above results to obtain rigorous and uniformly valid analyt- 
ical bounds on Nu{Ra) for convection with general Biot number thermal BCs. 

Theorem 4. For Rayleigh-Benard convection subject to thermal boundary conditions 
with Biot num,ber rj e [0, 00], let the control parameter R satisfy 

32 1 128 77 1 

i?>y for 77<-, i?>_^ for 77 > -. 

Then the Nusselt number is bounded according to 

Nu<^Ra^'\ (C5) 

independently ofrj. 
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Proof. We fix a Biot number 77 e [0, 00] and control parameter i? > 0, and define 
(5s(77, R), ba{r}, R) and i?c = 6s-R/ (6s — 1) as in LemmalU the constraint on R ensures tliat 
(5s(77, R) < 1/2. Then using S = Sg, we apply the background method with the piecewise 
linear background ts^, which is admissible by Lemma |31[b), so that by p.l5P " (|3.17p . 
A/'pwi,j)(^s, ^s) and I?pwi,j)(<Js, ^s) are rigorous bounds (upper and lower, respectively) for Nu 
and AT. From this theorem it also follows that 6s coincides with the balance parameter 
6c (|C 63(77, i?) = bc{r] , 6s{r] , R)) , so that we may apply the results of Lemma [U and 
conclude that Nu < A/'pwi,,,(<5s, 63) < 3/2Ss and AT > I?pwi, 77(^3, 6s) > (6s — l)/2, from 
which we may deduce a lower bound on Ra = RAT. 

In evaluating (|C 3|) at ^ = ^s, we consider the three cases in the definition of 6s('7, R) 
and ds{r], R): 

I: For < ?7 < 1/2 and R < Rsirj) = (8/3)rr^ in which case 6s - 1 = 1/2 and 
Ss = 2(2/3)i/2i?-i/2 > r,: 

ATy < A ^ Ra ^ RAT > ]r, 

and so Nu < |V6(4i?a)i/2 ^ ^^Ra^/^, 

II: For < 77 < 1/2 and R > Rs{r]), for which 6s - 1 = 5s/2f] and ^s = 2(r;/3)i/3^"i/3: 

iV^i < 4- = - f- V^' = ^AT >^R=l(^] i?2/3, 

- 2^s 4 V?7/ " 4r/ 2 V3r/2 ^ 

which gives R^^^ < 21/231/6,^1/3^^1/2 ^^j^^ ^^lus Nu < f^/6 iJa^/^ 

IIL For 1/2 < ?7 < 00 and i? > 0, for which 6s - 1 = (1 + 2ri)5s/^ri and 5s = 
2(2/3)i/3(,;/(l + 277))i/3i?-i/3: 

- 2<5s 4 V2 77 ; \ 7^ J 



and so < 25/631/6(77/(1 + 2?7)) 1/3 iZa^/^ and A^m < |\/6 iia^/^. 



□ 



Remark: The upper bound (|C 5P A^u < Crig,unif^a on Nu{Ra), valid uniformly in i? 
(and hence i?a) and in 77 g [0, cx)], was obtained at the cost of weakening the prefactor 
Crig,unif — 3\/6/4 « 1.837 relative to the asymptotic scaling (|5.15p Nu < Casym^a^^^, 
C'asym = \/2/27 ~ 0.272 for 7/ > 0. As in Lemma [2l the prefactor C — C^g in the rigorous 
bound Nu < CRa^^^ may be improved for particular (intervals of) rj. For instance, for 
fixed temperature BCs (77 = 0), (|5.3p establishes the bound with C ~ 3\/6/32 « 0.230; 
for fixed flux BCs (77 = 00) one can prove that it is sufficient to take C — •\/5/54 k, 0.304; 
and uniformly for 1/2 < 7; < 00, the use of the estimates (IC 41) from Lemma [5] instead of 
(|C 3p immediately allows one to improve the prefactor in (|C Sp to C = 5\/2/12 ~ 0.589. 
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